I -  :  l 

^AD-A0®4  345 
UNCLASSIFIED 

OHIO  STATE  UNIV  COLUMBUS  ELECTROSCIENCE  LAB  F/6  20/14 

RADIATION  BY  SOURCES  ON  PERFECTLY  CONDUCTING  CONVEX  CYLlNptRS  «— ETCCU) 

JAN  BO  L  ERSOY*  P  H  PATHAK  F19628-77-C-0107 

ESL-784641-2  RADC  -TR-79-2B1  NL 

n 

H 

. 

i 

j 

1 

■ 

_  f 

mm 

■■■ 

■■ 

■IB 

mm 

mm1 

mm 

m^m 

mm 

■ _ 

vs}  Cf'tU  ADA084345 


RADIATION  BY  SOURCES  ON 
PERFECTLY  CONDUCTING  CONVEX 
CYLINDERS  WITH  AN  IMPEDANCE 
SURFACE  PATCH 

The  Ohio  Stolo  University  ElectroScience  Laboratory 


L.  Ersoy 
P.  H.  Pathak 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED 


ROME  AIR  DEVELOPMENT  CENTER 

Air  Force  Systems  Command 

Griffiss  Air  Force  Base,  New  York  13441 


I 


80  5  X®  0^1 


I 


This  reporc  has  been  reviewed  by  the  RADC  Public  Affairs  Office  (PA) 
and  is  releasable  to  the  National  Technical  Information  Service  (NTIS) .  At 
NTIS  it  will  be  releasable  to  the  general  public,  including  foreign  nations. 

RADC-TR-7 9-281  has  been  reviewed  and  is  approved  for  publication. 


APPROVED: 


NICHOLAS  P.  KERNWEIS 
Contract  Monitor 


APPROVED: 


ALLAN  C.  SCHELL,  Chief 
Electromagnetic  Sciences  Division 


FOR  :OMMANDEr 


"JOHN  P.  HUSS 

Acting  Chief,  Plans  Office 


If  your  address  has  changed  or  if  'u 
mailing  list,  or  it  the  addressee 
please  notify  RADC  EEA)  Hanscom  A”  ' 
maintaining  a  curr  at  failing  list. 


wish  to  be  removed  from  the  RADC 
">  longer  employed  fcy  your  organisation 
A  01731  This  will  assist  us  In 


Do  not  return  this  copy.  Rete’.n  or  destroy. 


FICATION  OF  THIS  RAGE  (Whtn  Dmtm  Bnttrtd) 


EPORT  DOCUMENTATION  PAGE 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


RECIPIENT'S  catalog  number 


FORMING  ORGANIZATION  NAME  ANO  ADDRESS 


Ohio  State  University  ElectroScience  Laboratory 
Department  of  Electrical  Engineering 
Columbus  OH  43212  f  lL  l 


<1.  CONTROLLING  OFFICE  NAME  ANO  ADDRESS 


Deputy  for  Electronic  Technology  (RADC/EEA) 
Hanscom  AFB  MA  01731 


I  MONITORING  AGENCY  NAME  A  ADDRESS/!/  dltl»rwnt  from  Controlling  Oliicm)  IS.  SECURITY  CLASS,  (oi  thla  tmporl) 


sinzfl 


UNCLASSIFIED 


IS*.  DECLASSI  FI  CATION/ DOWN  GRADING 
SCHEDULE 

N/A 


is.  DISTRIBUTION  STATEMENT  ( a I  thlm  Rm, 


Approved  for  public  release;  distribution  unlimited. 


rfi/zr* 


IT.  DISTRIBUTION  STATEMENT  /«/  (A*  mbmlrmc t 


is  supplementary  notes 


1 


The  material  contained  in  this  report  is  also  used  as  a  dissertation  submitted 
to  the  Ohio  State  University  Graduate  School  as  partial  fulfillment  for  the 
degree  Doctor  of  Philosophy. 

RADC  Project  Engineer:  Nicholas  P.  Kernweis  (EEA) 


19.  KEY  WORDS  (Continue  on  rovoroo  otdo  it  necoooary  end  Identify  by  block  numbor) 


Impedance  Surfaces 
Perfectly  conducting  cylinders 
partly  covered  by  an  Impedance  patch 
Radiation  from  Sources  on  an 
impedance  patch 


Asymptotic  high  frequency  solutions 
Currents  excited  by  sources  on  curved 
Impedance  surfaces 

Perturbation  type  integral  equations 
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^This  report  deals  with  an  asymptotic  high  frequency  analysis  of  the  radiation 
patterns  of  a  magnetic  line  source,  or  a  magnetic  line  dipole,  located  on  a 
uniform  impedance  surface  patch  which  partly  covers  an  electrically  large 
perfectly  conducting  convex  cylinder.  This  work  is  relevant,  for  example,  to 
the  analysis  of  fuselage  mounted  airborne  antennas  for  satellite  communication 
purposes.  In  the  latter  application  impedance  surface  patches  may  be  employed 
to  increase  the  radiation  intensity  near  the  horizon  over  that  which  would 
ordinarily  exist  in  the  absence  of  any  such  impedance  loading.  In  the  present 
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^analysis,  the  impedance  surface  patch  Is  represented  as  an  equivalent  aperture 
in  the  rest  of  the  perfectly  conducting  convex  cylinder.  Approximate 
asymptotic  expressions  for  the  "equivalent  aperture  distribution;"  i.e.,  for 
the  "currents  on  the  impedance  surface  patch"  are  developed  in  this  work  for 
the  two  source  types;  these  expressions  for  the  currents  are  valid  in  the 
neighborhood  of  the  source.  The  radiation  pattern  of  this  configuration  is 
then  .found  in  a  straightforward  manner  by  numerically  Integrating  this  current 
distribution  in  conjunction  with  a  simple  and  accurate  asymptotic  high 
frequency  form  of  the  perfectly  conducting,  convex  cylinder  Green's  f unction. _ 
Since  the  cylinder  Green's  function  is  employed,  one  needs  to  integrate  only 
over  the  currents  which  exist  on  the  Impedance  surface  patch  (i.e.,  over  the 
equivalent  aperture)  and  not  over  the  currents  on  the  remaining,  perfectly  ,/ 
conducting,  portion  of  the  cylinder .J^Radiat ion  pattern  calculations  based  on 
this  analysis  are  found  to  compare  quite  well  with  a  numerical  moment  method 
type  solution,  a  previous  GTD  solution  for  those  impedances  which  support  an 
Elliot  mode,  and  also  with  some  presently  available  experimental  results. _ It 
is  seen  that  a  moderately  large  increase  in  the  radiation  intensity  can  be'" — 
obtained  in  the  vicinity  of  the  horizon  for  cylinders  loaded  by  an  impedance 
surface  which  support  Elliot  type  modes. 
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EVALUATION 

This  Is  the  final  report  on  Contract  FI 9628-77 -C-01 07 .  It  covers 
research  and  experiments  performed  during  the  period  6  January  1977  to 
5  January  1979.  The  objective  of  the  work  was  to  Increase  the  radiation 
Intensity  near  the  horizon  of  a  magnetic  line  source,  or  a  magnetic  line 
dipole,  located  on  a  uniform  Impedance  surface  patch  which  partly  covers 
an  electrically  large  perfectly  conducting  convex  cylinder.  Various 
methods  for  analyzing  and  computing  the  radiation  patterns  of  such  magnetic 
sources  are  described  and  the  near  and  far  fields  can  be  predicted  fairly 
well . 
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CHAPTER  I 
INTRODUCTION 


An  approximate,  asymptotic  high  frequency  analysis  of  the 
radiation  patterns  is  presented  for  a  two-dimensional  (2-D)  mag¬ 
netic  line  source  or  a  magnetic  line  dipole  source  located  on  a 
uniform  Impedance  surface  patch  which  partly  covers  an  electrically 
large,  perfectly  conducting  convex  cylinder  as  in  Figure  1.  The 
purpose  of  this  work  Is  to  assess  the  effects  of  such  an  impedance 
loading  on  the  radiation  patterns  with  a  view  towards  being  able 
to  Increase  the  pattern  level  In  the  vicinity  of  the  horizon. 

The  latter  application  is  of  interest  in  the  design  of  a  crossed- 
slot  type  fuselage  mounted  airborne  antenna  for  satellite  communi¬ 
cation  purposes.  On  a  large,  perfectly  conducting  convex  cylinder, 
a  circumferential  slot  typically  radiates  a  field  at  the  horizon 
(or  shadow  boundary)  which  is  roughly  20  dB  below  that  of  an  axial 
slot.  In  the  crossed  slot  application,  it  is  therefore  desirable 
to  enhance  the  fields  of  a  circumferential  slot  near  the  horizon 
by  impedance  loading,  while  at  the  same  time  not  significantly 
degrading  the  performance  of  the  axial  slot  in  the  crossed  slot 
antenna  configuration. 

In  the  present  method  of  analysis,  the  impedance  surface  is 
represented  by  an  equivalent  aperture  in  an  otherwise  perfectly 
conducting  convex  cylinder.  Approximate,  asymptotic  expressions 
for  the  "equivalent-aperture  distribution"  corresponding  to  the 
surface  currents  excited  by  a  source  on  the  impedance  patch  are 
developed  for  both  the  magnetic  line  source  and  the  magnetic  line 
dipole  cases.  A  magnetic  line  source  simulates  a  thin  axial  slot; 
whereas,  a  magnetic  line  dipole  simulates  a  thin  circumferential 
type  slot.  The  asymptotic  expressions  obtained  for  the  surface 
currents  on  the  impedance  patch  remain  valid  in  the  vicinity  of 
the  source  at  Q'.  The  radiation  pattern  of  the  configuration  in 
Figure  1  is  then  found  in  a  straightforward  fashion  by  numerically 
integrating  this  current  distribution  in  conjunction  with  an  ac¬ 
curate,  asymptotic  form  of  the  perfectly  conducting,  convex  cyl¬ 
inder  Green's  function  as  given  by  Pathak  and  Kouyoumjianfl] ;  a 
slightly  improved  version  of  this  Green's  function  also  appears 
in  a  report  by  Pathak  and  Huang[2j.  As  a  result  of  employing  this 
special  cylinder  Green's  function  rather  than  the  usual  free  space 
Green's  function,  one  needs  to  integrate  the  currents  only  over 
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the  Impedance  surface  patch.  This  method  of  analysis  will  be  re¬ 
ferred  to  as  the  surface  current  Integration  method  (SCIM).  Radi¬ 
ation  pattern  calculations  based  on  this  SCIM  are  found  to  compare 
quite  well  with  those  obtained  via  a  numerical  moment  method  type 
solution,  a  previous  GTD  solution  for  those  Impedances  which  sup¬ 
port  an  Elliot  mode, and  also  with  some  presently  available  experi¬ 
mental  results.  A  moderately  large  Increase  In  the  radiation 
Intensity  can  be  obtained  in  the  vicinity  of  the  horizon  for  cyl¬ 
inders  loaded  by  an  Impedance  surface  which  support  Elliot  type 
modes[3].  An  Elliot  mode  field  Is  one  which  propagates  on  a  curved 
Impedance  boundary  with  negligible  leakage,  and  Is  therefore  usually 
the  dominant  contributor  to  the  field  on  such  a  surface;  further¬ 
more,  as  the  curvature  of  the  surface  vanishes,  it  reduces  to  the 
usual  "bound"  surface  wave  mode  on  a  planar  impedance  surface. 

It  Is  noted  that  an  Elliot  mode  can  be  excited  by  a  magnetic  line 
source  If  the  surface  impedance  Is  inductive;  whereas,  it  can  be 
excited  by  a  magnetic  line  dipole  source  If  the  surface  Impedance 
Is  capacitive. 


IMPEDANCE  SOURCE 

PATCH  /  OBSERVATION  POINT 

^  ON  IMPEDANCE  PATCH 


Q'Q,  *  L,  (ARC  LENGTH) 
0^*  L2  (ARC  LENGTH) 
L  (LENGTH  OF  IMPEDANCE 
PATCH)  *  L,  +  Lt 


Figure  1— An  Illustration  of  an  antenna  on  a  perfectly  conducting 
convex  cylinder  which  is  partly  covered 
by  an  Impedance  patch. 

This  present  work  represents  an  Important  extension  and  gen¬ 
eralization  of  some  of  the  earlier  work  02,4,5],  The  earlier  ap¬ 
proach  by  Pathak  [4  ]  which  led  to  a  preliminary,  asymptotic  solution, 
for  this  problem  could  be  designated  as  the  diffraction  coefficient 
or  the  geometrical  theory  of  diffraction  approach  (GTD).  A  second 
approach  presented  by  Pathak  and  Huang  D?]  which  is  based  on  the 
numerical  solution  of  the  Integral  equation  for  this  problem  is 
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designated  as  the  MM-GTD  approach  [2] .  This  numerical  solution 
was  based  on  the  moment  method  with  some  important  simplifications 
resulting  from  the  use  of  uniform  and  accurate  GTD  approximations 
for  the  excitation  and  the  kernel  of  the  integral  equation.  This 
MM-GTD  solution[2]  is  employed  primarily  as  an  independent  check 
on  the  GTD  solution  ft];  it  also  serves  as  a  check  on  the  SCIM  solu¬ 
tion.  The  MM-GTD  solution  in  [21,  which  employs  the  cylinder  sur¬ 
face  Green's  function  as  a  kernel,  is  rather  well  suited  for  the 
TE  case.  However,  one  needs  to  be  careful  in  employing  it  to 
the  TM  case  because  the  kernel  is  far  more  singular  than  for  the 
TE  case,  and  it  then  appears  to  be  necessary  to  properly  incorpo¬ 
rate  the  edge  conditions  on  the  fields  at  Q,  and  in  order  to 
obtain  a  completely  stable  numerical  result.  While  some  initial 
success  has  been  obtained  for  the  TM  case,  further  work  is  neces¬ 
sary  to  complete  the  MM-GTD  treatment  more  satisfactorily  for  this 
case.  Preliminary  modifications  and  extensions  of  this  work  in 
[2]  and  [4]  were  reported  later  by  Pathak  and  Ersoy  [5]  in  which 
the  SCIM  procedure  was  outlined  and  a  study  of  the  asymptotic  ap- 
roximations  for  the  surface  currents  to  be  employed  in  this  SCIM 
procedure  was  initiated. 

The  work  based  on  the  GTD  approach  in  [4]  is  an  extension 
of  the  GTD  solution  of  Hwang,  Kouyoumjian  and  Pathak [6]  for  anal¬ 
yzing  the  radiation  by  a  magnetic  line  source  on  a  perfectly  con¬ 
ducting  box  that  is  partly  covered  by  a  finite  dielectric  panel 
(over  the  face  on  which  the  source  is  located).  In  contrast  to 
the  work  of  Hwang  et  al.[6],  a  GTD  analysis  of  the  curved  geometry 
in  Figure  1  requires  the  solutions  to  more  and  different  canonical 
problems,  because  of  the  complexity  introduced  into  the  analysis 
by  the  presence  of  the  curvature.  In  this  context,  it  is  noted 
that  the  GTD  solutions  including  the  one  in  [1]  and  [4]  are  based 
on  the  local  properties  of  wave  propagation  and  diffraction  at 
high  frequencies,  and  can  hence  be  "built  up"  from  the  asymptotic 
solutions  to  appropriate  canonical  problems  which  accurately  model 
the  local  regions  from  where  excitations  or  diffractions  of  the 
fields  take  place.  For  example,  the  field  at  a  far  zone  observa¬ 
tion  point  which  is  generated  by  the  source  at  Q'  is  obtained  by 
summing  up  the  fields  which  propagate  along  the  direct  ray  path, 
from  the  source  to  the  field  point  as  well  as  along  the  various 
diffracted  ray  paths  such  as  those  emanating  from  the  ends  of  the 
patches  at  Q.  and  Q?,  and  also  those  that  are  shed  tangentially 
after  creeping  around  the  back  of  the  cylinder  upon  being  launched 
from  Q.  and  Q-.  The  ray  picture  of  the  GTD  as  illustrated  in  Figure 
2  is  pnysicany  appealing  as  it  provides  an  insight  into  the  nature 
of  the  radiation  and  diffraction  mechanisms.  The  GTD  analysis[4] 
was  restricted  to  an  impedance  surface  which  would  support  an  Elliot 
mode  and  for  which  this  mode  alone  would  be  the  dominant  contributor 
to  the  fields  diffracted  from  the  edges  of  the  impedance  patch 
at  Q,  and  Q«.  The  other  component  of  the  surface  field  (i.e., 
besides  the  Elliot  mode  field)  diffracts  not  only  from  the  edges 
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2— Rays  and  shadow  boundaries  associated  with  the 
GTD  solution  of  [4]. 
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of  the  impedance  patch  but  also  from  the  impedance  surface  via 
a  continual  leakage  along  the  ray  propagation  path  from  Q*  to  Q, 
or  Q2.  If  this  leakage  is  significant,  then  the  diffraction  from 
just^the  Elliot  mode  may  not  provide  a  sufficiently  accurate  radi¬ 
ation  pattern  near  the  horizon.  One  expects  this  leakage  effect, 
which  is  partly  governed  by  the  type  of  impedance,  to  become  im¬ 
portant  as  the  cylinder  size  decreases.  Furthermore,  if  the  type 
of  impedance  chosen  does  not  support  an  Elliot  mode  field,  then 
this  other  field  component  is  all  that  is  left.  For  the  latter 
case,  the  GTD  analysis  of  UJ  must  be  substantially  modified  to 
obtain  a  uniform  diffraction  coefficient  for  the  diffraction  from 
the  edges  of  the  patch,  of  this  component  of  the  field  which  con¬ 
tinually  leaks  (or  sheds)  energy  from  the  impedance  surface.  A 
uniform  asymptotic  solution  of  an  appropriate  canonical  problem 
which  would  lead  to  such  a  uniform  diffraction  coefficient  is  some¬ 
what  complicated  and  is  not  attempted  at  the  present  time.  The 
GTD  technique  presented  in  [4]  is  thus  useful  for  sufficiently 
large  cylinders  with  a  sufficiently  long  impedance  patch  on  which 
an  Elliot  mode  can  establish  itself.  It  is  noted  that  a  pure  Elliot 
mode  is  established  beyond  a  certain  minimum  launching  distance  on 
the  impedance  patch.  For  more  general  situations;  namely,  for 
the  case  when  the  cylinders  are  only  moderately  large  and  for  im¬ 
pedance  surfaces  which  are  not  necessarily  restricted  to  those 
that  can  support  an  Elliot  mode  field,  the  approach  based  on  the 
SCIM  appears  to  lead  to  a  relatively  simpler  and  sufficiently  prac¬ 
tical  method  of  solution  for  the  problem  in  Figure  1.  Furthermore, 
this  SCIM  solution  is  constructed  without  essentially  losing  much 
of  the  physical  insight  provided  by  GTD  solutions.  Also,  the  SCIM 
can  take  into  account  impedance  surface  patches  which  are  reasonably 
small;  more  will  be  said  about  this  in  Chapter  III.  It  is  noted 
that  while  the  SCIM  can  handle  more  general  situations  than  the 
previous  GTD  solution,  it  is  also  restricted  at  this  time  as  is 
the  latter  solution  to  an  impedance  surface  which  is  uniform  (or 
at  the  most  one  which  is  very  slowly  varying  if  it  is  non-uniform). 
The  MM-GTD  procedure[2]  can  taki  non-uniform  impedance  into  con¬ 
sideration  quite  easily;  however,  it  does  not  provide  the  physical 
insight  afforded  by  the  SCIM. 

It  may  >e  mentioned  that  Wait[7-12]  has  done  considerable 
work  on  the  asymptotic  calculation  of  the  surface  fields  which 
are  diffracted  past  an  impedance  discontinuity  as  at  Q,  and  Qp 
of  Figure  1.  Indeed,  some  of  his  work[7,9,12]  has  been  helpful 
in  the  asymptotic  analysis  of  the  currents  which  exist  on  a  curved 
impedance  boundary.  Nevertheless,  Wait's  results  have  not  been 
directly  applied  here,  as  they  do  not  furnish  surface  current  repre¬ 
sentations  which  are  suitable  in  the  construction  of  an  approximate 
solution  based  in  particular  on  the  SCIM.  Additional  comments 
comparing  Wait's  results  with  the  ones  developed  here  are  presented 
in  Section  I  IB .  Also,  Shapira,  Felsen  and  Hessel[13,14]  have  anal¬ 
yzed  the  radiation  from  an  active  slot  in  an  array  of  slots  on 


a  perfectly  conducting  convex  cylinder  via  asymptotic  techniques; 
they  employ  a  local  periodic  structure  (LPS)  model  in  which  the 
array  is  viewed  as  an  effective  surface  impedance.  They  also  in¬ 
clude  the  edge  effects  arising  from  finiteness  of  the  array.  It 
is  noted  that  the  surface  fields  obtained  in  the  LPS  method  of 
[13]  is  not  valid  in  the  close  neighborhood  of  the  source;  whereas, 
in  the  SCIM,  it  is  important  to  obtain  a  surface  field  represent¬ 
ation  which  remains  valid  in  the  source  region  as  is  done  in  this 
work.  Furthermore  the  radiated  field  in  the  LPS  method  of  [13]  is 
not  as  easily  obtained  in  the  vicinity  of  the  horizon  as  it  is 
with  the  SCIM. 

The  organization  of  this  dissertation  is  as  follows.  Chapter 
II  describes  the  method  of  solution  based  on  SCIM.  The  starting 
point  of  this  solution  is  a  generalized  reaction  formulation  of 
the  problem  which  leads  to  an  integral  representation  for  the  fields 
radiated  by  the  configuration  in  Figure  1.  This  integral  repre¬ 
sentation  is  obtained  in  terms  of  the  surface  fields  on  an  equiva¬ 
lent  aperture  corresponding  to  the  extent  of  the  impedance  surface 
patch.  The  approximate,  asymptotic  expressions  developed  in  here 
for  these  surface  fields  over  the  impedance  patch  are  also  indi¬ 
cated  in  that  chapter  together  with  a  brief  description  on  the 
final  steps  involving  the  integration  of  these  surface  fields  to 
obtain  the  radiated  fields.  Chapter  III  discusses  the  accuracy 
of  the  solution  obtained  via  SCIM.  Numerical  results  illustrating 
the  effect  of  different  values  of  impedances  and  cylinder  sizes 
on  the  radiation  patterns  of  a  magnetic  line  source  or  a  magnetic 
line  dipole  source  at  Q1  are  presented  in  Chapter  IV.  The  major 
conclusions  of  this  work  are  presented  in  Chapter  V.  Analytical 
details  pertaining  to  the  development  of  the  present  SCIM  solution 
are  outlined  in  several  appendices. 
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CHAPTER  II 
METHOD  OF  SOLUTION 


The  method  of  solution  employed  here  for  treating  the  radi¬ 
ation  problem  depicted  in  Figure  1  is  discussed  in  this  chapter. 

As  pointed  out  earlier,  in  this  solution  based  on  the  SCIM,  the 
impedance  surface  patch  is  viewed  as  an  equivalent  aperture  in 
the  rest  of  the  perfectly  conducting  convex  cylinder.  The  field 
radiated  by  this  configuration  is  then  expressed  in  terms  of  an 
integral  just  over  the  surface  current  distribution  which  exists 
on  the  aperture  (or  the  impedance  surface  patch)  because  one  em¬ 
ploys  an  asymptotic  high  frequency  form  of  the  perfectly  conducting 
convex  cylinder  Green's  function  in  this  radiation  integral.  The 
formulation  of  this  radiation  integral  over  the  equivalent  current 
distribution  in  the  aperture  is  discussed  in  part  A  of  this  chapter. 
The  surface  current  distribution  on  the  impedance  patch,  which 
is  produced  by  the  source  at  Q*  is  approximately  found  from  the 
asymptotic  solutions  to  two  auxiliary  canon-.al  problems.  The 
solutions  to  these  auxiliary  canonical  problems  and  the  form  of 
the  surface  currents  on  the  impedance  patch  are  discussed  in  part 
B.  The  final  step  in  the  solution  involves  incorporating  the  sur¬ 
face  currents  found  from  part  B  into  the  radiation  integral  of 
part  A,  and  then  evaluating  this  integral  to  obtain  the  radiated 
field.  The  latter  integration  is  performed  efficiently  and  in 
a  straightforward  manner  via  standard  numerical  techniques,  as 
discussed  in  part  C.  In  the  following  analysis,  an  eJa)t  time 
dependence  will  be  assumed  and  suppressed. 

A.  Formulation  of  the  Radiation  Integral 

The  pertinent  radiation  integral  may  be  obtained  in  several 
different  ways.  The  procedure  by  which  this  radiation  integral 
is  obtained  in  the  present  work  is  based  on  a  generalization  of 
the  reaction  theorem[l5,16j.  The  same  radiation  integral  could 
also  be  obtained  via  an  application  of  the  scaler  Green's  theorem 
(since  the  problem  considered  here  is  two-dimensional),  or  the 
dyadic  Green's  theorem  (and  it's  specialization  to  the  two-dimen¬ 
sional  vector  problem  considered  here),  or  an  application  of  the 
compensation  theorem  for  electromagnetic  fieldsLl7J.  In  the  methods 
based  on  the  Green's  theorems,  an  appropriate  Green's  function 
which  satisfies  the  boundary  conditions  for  a  perfectly  conducting 
convex  cylinder  must  be  employed.  The  generalized  reaction  theorem 
which  is  employed  here  for  formulating  the  radiation  Integral  is 
somewhat  related  to  the  other  methods  mentioned  above;  it  is  chosen 
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7x1^  =  -JkZ0T^;  vxl^  =  jkY0^  +  7t;  if  (la,b) 

and 

VX^  =  -Jkz0^  -  Hj,-  vxI^  =  jkY0^;  If  T-Bt  (2a, b) 

Furthermore,  (7  ,H )  satisfy  the  radiation  condition  over  the  surface 
Z  at  infinity,  and°E  satisfies  the  following  boundary  conditions 
on  the  surface  s;  0 

fix^’m  •  0  ,  on  s.  (3) 

In  (1)  and  (2)  k  is  the  free  space  wavenumber,  and  Z0=1/Yq 
is  the  free  space  Impedance  (i.e.,  ZQ=120 *  ohms). 

The  fields,  T.,  TL  of  the  perturbed  problem  also  satisfy  the 
Maxwell's  equations 


VxFb  =  -JkZ0Hb-M;  vxTTb  =  jkY^ 


(4a, b) 


and  the  radiation  condition  over  the  surface  z  at  infinity.  Also 
Fb  and  flb  satisfy  the  following  boundary  conditions  on  s: 

\  -  (n«fb)n  =  ZsnxTTb^|  on  Sj  (5a) 


nxrb  =  0 


s*sl+s2 


where  Z  denotes  the  value  of  the  surface  impedance  in  the  perturbed 
problem. 

Combination  of  (2),  (4)  and  (1),  (4)  respectively  yields 

»•  -  v1?]  *  'K-W  TVcase  <«> 

and 

7’^oxffb  '  =  ™z‘case  (7) 

Applying  the  divergence  theorem  to  (6)  and  (7)  and  using  the  radiation 
condition  on  I,  and  the  boundary  conditions  (3)  and  (5)  yields: 


,  -  ,  . . >r,  *  ir  * 

-  *  -  •*  '<  .  ■ 
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y  dsflJ(Q|P').[Zs(a)fixfixHb(0|P)] 


///  dvCHS(P|P* )*W<P)^b(P* |P)-T(P*) 
v  If  T=X 


///  dvDff(P|P,).M(P)-1L(P'|P).T(P')]. 
.  v  if  T=R. 


(8) 


The  volume  v  Is  enclosed  between  the  closed  surfaces  s*s.+s2  and  z 
(which  receeds  to  infinity).  The  notation  of  the  type  (A|Bf  In 
the  argument  of  the  field  functions  and  (T.  ,TT.  )  Is  Introduced 
for  convenience,  and  it  refers  to  the°f1eld  at  A  doe  to  source 
at  B.  Specializing  the  above  equations  to  the  two-dimensional 
field  and  source  configurations  (i.e.,  configurations  which  are 
uniform  in  the  2-direction),  one  obtains 


/  d^(Q|P').[Z  (Q)ftxfixWb(Q|P)]  = 


//dxdyR®(P|P' )  •H(P)+//dxdyr.(P,|P)  -T(P-) 

SM  ST 

M  T  if  T=3t, 


JJdxdy(^(P |P* )  .M(P)-//dxdyffb(P*  |P)  -T(P' ) 


SM  ’  ST 


if  T-Ht. 


(9) 


The  regions  sM  and  sT  are  the  cross-sectional  areas  of  the  2-D  source 

distribution  M  and  T,  respectively.  Also,  ds  In  (8)  is  given  as  ds 
=  dj,dz>0  with  dz>0;  hence,  the  2-D  arc  length  integral  on  the  boundary 
s,  (LHS  of  (9))  Is  chosen  to  proceed  in  the  direction  for  which  dfi,  in¬ 
creases  in  the  positive  sense.  Let  the  2-D  sources  W  and  T  be 


{^fiCp-ppi)*  if 

(10b) 

mt6(p-Pp.).  if  T»Rt 

Let  the  point  P'  In  (9)  receed  to  Infinity;  whereas,  the  point  P  Is 
moved  down  to  the  point  Q*  on  the  surface.  Thus, 

/  d^(Q|P').[Zs(Q)fixnxT?b(Q|Q')]  -  ^(Q*  |P* )^b(P'  |Q* )‘3t 
S1 

for  the  TMZ  case  with  m*TMd  and  3t*z  (11a) 


/  dJlflJ(Q|P-)*[2s(Q)^xnxHb(QjQ')]  =  flJ(Q' |P' )»m-Hb(P' |Q')*mt 
S1 


for  the  TE2  case  with  m=zM  and  mt=z.  (lib) 


h ( ? ' I Q ' )  i°  the  fM  problem  pertaining  to  (11a)  Is  the  electric 
“which  Is  producedzat  P'  (in  the  far  zone)  by  a  unit  strength  tan- 


The 
field1 

gentlal  magnetic  line  dipole  source  at  Q'  on  the  surface  s.;  whereas, 
H.(P'|Q')  Is  the  magnetic  field  in  the  TE  problem  In  (lib)  which 
Is  produced  at  P'  (In  far  zone)  bv  a  un1tzstrength  magnetic  line 
source  at  Q‘  on  the  surface  s,.  H~  (Q*  or  QlP'j  and 
Q | P ' )  are  the  magnetic  f lei ds* induced  at  Q'  6r 
surface  s  of  Fltjyr'  “  J  *  "  "  - ^  " 


(In  far  zone) 

- "  (O'  or 

on  the  unperturbed 
jre  3  due  to  a  test  source  T  at  P*  In  the  far  zone, 
the  test  source  Is  an  electric  line  source; 


?|jemtest  source  is  a  magnetic  line  source.  The 


In  the  case  ol 
wheras,  for  Tr  , 

expressions  for  TP’"'  are  available  In  [1,2]  in  terms  of  the 
6TD  approx1mat1on°wh1ch  will  be  employed  here.  The  GTD  approxi¬ 
mation  Is  valid  for  electrically  large  cylinders;  the  latter  Is 
assumed  to  be  valid  In  the  present  analysis.  Thus,  the  F.  and 
H.  on  RHS  of  (11a)  and  (lib)  may  be  found  if  nxH.  (Q I Q* )  iRsidg 
tne  Integrals  on  the  LHS  of  (lla;b)  are  known.  The  field  nxfl.(Q|Q') 
Is  the  tangential  magnetic  field  at  Q  on  the  surface  s,  of  Figure 
3  due  to  a  source  at  Q*  also  on  s,.  In  the  TM  case,  the  source 
at  Q'  Is  a  tangential  magnetic  line  dipole;  whereas,  In  the  TE 
case,  the  source  Is  a  magnetic  line  source  as  mentioned  earlier. 


An  alternative  expression  for  F.(P'|Q')  for  the  TM  -case  of 
(11a)  can  be  expressed  In  terms  of  electric  field  quantities  alone 
via 
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(12) 


T^(Q‘|P*)  •  n  -  -IJ(P'IQ')  •  3t  with  3t  *  z 

as  follows: 

Ib(p,|Q')  *  2  =^(p,|Q')  *  z  +  jr  /  diC^(P-IQ)  •  z] 

d  S1 

•  Clb(Q|Q')  •  z  3  (13) 

Furthermore,  it  is  shown  below  in  section  B  that  the  unknown  I.  (Q|Q' ) 
can  be  expressed  as 

tb(Q|Q')  •  z  *  -Md  6(Q|Q')  +  Ib(Q+|Q')  •  z  (14) 

In  (14),  r.(QlQ')  •  z  is  the  tgngentjaj^  electric  field  of  the  mag¬ 
netic  11neDdipole  source,  *  t  Mdfi(p-p0,)  and  it  is  observed  on 
the  surface  s,  wh1ch„j1es  "Just  underneath"  the  magnetic  line  dipole 
source;  whereas,  T.  (Q  IQ')*z  Is  the  tangential  electric  field  of 
the  same  source,  and  tnls  field  is  observed  on  a  surface  sj  which 

lies  "just  above"  $£e  magnetic  line  dipole  source.  Here  (T  must 
be  Interpreted  as  Q  *  1 1m(p*p,-»-e;4>)  where  Q'  *  (p ' ,4> * )  and  e  Is 

£■*■0 

a  pgsltlve  number:,  Explicit  expressions  for  the  electric  field 
Fb(Q|Q')  and  Fb(Q  |Q‘)  of  (14)  are  given  in  Appendix  B.  Equation 
(14)  Is  a  well  known  expression  for  the  discontinuity  of  the  tan¬ 
gential  ejectrlc  field  as  one  moves  from  the  surface  s.  to  the 
surface  Sj  across  a  magnetic  current.  1 

When  (14)  Is  Incorporated  into  (13),  a  further  simplification 
results  so  that  the  final  form  of  ( P ’ | Q ' )  becomes 

Tb(p,lQ,)ar  /  d*Op‘!Q)  -  zJCrb(Q+|Q')  •  £]  (15) 

d  + 


Ijjj  (15)  It  Is  noted  that  the  Integral  Is  performed  on  the  surface 
Sj  "just  above"  the  line  dipole  source. 

The  Importance  of  Introducing  (15)  will  become  evident  In 
section  (Cohere  It  Is  pointed  out  that  the  expression  for  the 
unknown  e.  (Q  |Q')  In  the  TM  case  Is  less  complicated  than  the 
one  for  tRe  unknown  W^Q1  |Q;  of  (11a)  for  the  same  TMZ  case. 

An  approximation  for  the  unknown  Z  nxnxTT.  (Q(Q')  in  (llj;llb) 
will  be  Introduced  next  together  with  an  approximation  for  z  • 

rb(Q+|Q'). 
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Approximation  for  the  Surface  Currents 
on  the  Impedance  Patch 


The  current  Z 
length  Impedance 


Z  nxnxTT- (Q JQ • ) 
sarf ace  patch  s 


of  (11)  which  exists  on  the  finite 
j  may  be  expressed  as 


ZsnxnxHb(Q|Q') 


case 

case 


(16a) 


and  similarly  the  current  "E".  ( Q+ 1 Q * )  •  z  (=  x»+nxIb(Q+|Q* ))  of  (15) 
which  exists  on  the  impedance  surface  patch  s^  is, 

tb(Q+lQ')  •  z  =  -asMd[Ge(Q+|Q')  +  G®(Q+|Q* )]for  the  TM2-case 

(16b) 


where  a  is  defined  in  (20c).  M  or  M  .  are  the  strengths  of  the 
magnetic  line  source  or  line  dipole  source  at  Q'  for  the  TE  or 

§ 

TMz  case,  respectively.  The  - jkYQZs  1  1  ^ 


ource  ai 

{£) 


G  refers  to  that 


"d  J 

component  of  the  current  which  would  exist  if  the  cylinder  was 
completely  covered  by  an  impedance  surface.  The  quantity 

Gn  1 


-jkYZ 

os 


t-Mdl) 


to  the  current  -jkY  Z 
os 


s  then  viewed  as  the  correction  or  perturbation 

b 

G  which  arises  from  the 


truncation  of  the  impedance  surface  at  Q,  and  Q9.  The  term 
f  '  h 

G^  thus  refers  to  the  component  of  the  current  which 


-jkY  Z 
os 


WJ 


is  reflected  from  the  ends  of  the  patch  at  Qj  and  Q^. 


Explicit  expressions  for  Ge(Q+|Q')  is  given  in  (B-44a)  of 
Appendix  B.  A  comparison  Qf~{B7)  with  (B-44a)  reveals  that,  an 
asymptotic  expansion  for  Ge(Q+|Q')  can  be  obtained  from  the  asymp¬ 
totic  expansion  for  Gh  by  simply  changing  a.  and  B.  in  the  latter 
to  <*  and  Bc.  The  final  form  of  Gn  is  given  in  (17a) .  The  quantity 

o  "*x  *  u 

Gn (Q  |Q* )  which  is  analogous  to  G»  arises  from  the  fact  that  the 
impedance  patch  is  of  finite  extent,  and  therefore  there  are  re¬ 
flections  from  the  ends  of  the  patch  at  Q.  and  Q~. 
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In  particular  the  quantity  -jkY.Z,  J  Tg^  may  be  expllc- 

1'MJ 

Itly  found  from  the  canonical  problem  of  determining  the  surface 
fields  of  a  magnetic  line  source  or  a  magnetic  line  dipole  on  a 
circular  cylinder  which  Is  completely  covered  with  an  Impedance 
surface  as  In  Figure  4.  The  radius  of  the  cylinder  Is  denoted 
by  a.  The  points  on  the  Impedance  patch  are  specified  by  the  arc 
length  t(=a|$|)  along  the  surface  of  the  cylinder  and  It  Is  meas- 
ured.from  the  source.  In  appendix  B,  an  asymptotic  approximation 

to  Gs  for  the  circular  cylinder  case  Is  presented;  this  asymptotic 
approximation  Is  valid  for  electrically  large  cylinders  which  are 
of  interest  in  this  study.  While  an  exact  circumferential  eigen¬ 
function  (or  modal)  series  expansion  exists  for  G5  pertaining  to 
the  circular  cylinder  case,  It  is  poorly  convergent  for  large  cyl¬ 
inders.  Furthermore,  this  exact  model  series  representation  for 
the  circular  case  cannot  be  generalized  to  an  arbitrary  convex 
cylinder..  In  contrast,  an  asymptotic  high  frequency  approximation 

for  the  Gs  pertaining  to  a  circular  cylinder  becomes  Increasingly 
accurate  with  an  increase  In  the  cylinder  radius.  Also,  this  as¬ 
ymptotic  solution  can  be  generalized  rather  easily  via  the  local 
properties  of  high  frequency  field  solutions  to  treat  an  arbitrary 
convex  cylinder  with  a  slowly  varying  curvature  as  In  the  geometri¬ 
cal  theory  of  diffraction  (GTO)  procedure.  The  asymptotic  approxl- 

*5  i 


matlon  to  G3  for  the  circular  case  is  developed  in  sufficient  detail 
In  Appendix  B;  hence,  only  the  final  results  are  summarized  below. 


Gh^  1 


—0 

-bfkt 


h  H<2>(kt)F;<C)  -J  £  e'J6ht  eh  -  ^  e--"11  q( 


-*|f. 

4a  I  L 


,-jkt  J 


1+Jkt 


-Joo  5 


-jkt  »  --p 

'2irJ  ^ji  _I 


P=0  crp+(ahk”1m)2 


kfif 

\ 

+  1 


?  + ’}♦  k) F*(kt)]  b  ■ 


t>*st 


m 


t<%tm  (17a) 


(17b) 


and 
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_ 2  -  DIRECTED 

JT  MAGNETIC  LINE 
SOURCE 


(a)  te2-case 


«/>'-  DIRECTED 
MAGNETIC  LINE 
DIPOLE  SOURCE 


Q 


(b)  TMZ  -CASE 


e  4— Geometries  of  canonical  problems  of  a  line  source  on 
a  circular  cylinder  covered  completely  by  an 
impedance  surface. 
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,  Hf2)(kt)  a 

7J—R - F1(S)_J  7 


‘s*  . 


4-e  se+-s6(t)e 

kz0$  8s  s 


"hi1 


d(kt)' 


.[H(2)(kt){l-2z  ez2Q(z)}]  + 


t<»stm  (18a) 


2  -jop£ 

-*.1  V  .'**  f  «  .!  ;  ;  t>Wm 

4irk2  m  P=0  ap-(ask  m) 


(18b) 


where  t  is  the  maximum  total  impedance  patch  size  beyond  which 
the  asymptotic  representation  of  the  currents  in  (17a)  and  (18a) 
cannot  be  trusted.  An  empirical  equation  for  t  is  presented  in 
Chapter  III.  Also,  in  (17) 

F}<t>-l-<?  E9/2*-  (Wa) 


Fj(kt)  = 


1 

t  =  a)4»|  ;  £=  m  |$|  ;  m  =  (^) 


0,1  =  "Jk  Zo  5  8h  =^k2+ah  =  2irXb 


K  srJi;  ah =  'l«h|; 


-  - 

,  l  ,  ^  A 

=  o  *  ah<®  * 


(19b) 


(19c, d,e) 


(19f,g) 


(19h,i) 


( 19  j ) 


ap  are  the  roots  of  wi^-oLmk"  Wp(o)  =  0  where  Wp(o)  =  ir[Bi(o)-jAi(o)] 
is  the  Fock-type  Airy^f unction.  Tne  prime  over  Wp^indicates  the  deriv¬ 
ative  with  respect  to  the  argument.  These  roots  nave  been  extensively 
studied  by  Wait [7]  and  Logan [20],  (also  see  Appendix  B). 

Similarly,  the  quantities  in  (18}  are  defined  as  follows: 


Ff(C)  - 

i 

£/2 

5s  ♦  FS  ♦ ... 

(20a) 

-jest 

-jkt 

_2 

F|(kt)= 

e 

*s 

es  +  1 

3 

k  Jv 

Q(z)ez  ;  a$<0 

(20b) 

Ys 

°s  =  Y~  * 

0 

-y? 

>  2 
‘+as  = 

=  2*XS 

(20c, d) 

K  “S 

bs  =  IT 

U: 

“s 

-I«sl 

;  z  =  jbjkt 

(20e,f,g) 

in  which  the  z  of  (20-g)  is  not  to  be  confused  with  the  z  coordinates 
(the  present  analysis  is  independent  of  the  z  coordinate)  and  Y  =_1/Z 
is  the  surface  adnittance  in  mhos.  o„  are  the  roots  of  wl(o)  -of.k"im  s 
w2(a)  =0.  p  Z  s 

The  function  Q  is  related  to  the  complementary  error  function, 
erfc  such  that 

r—  00  2 

Q(R)  =  erfc(R)  =  /  e'x  dx  . 
c  R 

It  is  important  to  point  out  that  the  asymptotic  expansion  for 
G*  given  in  (18a)  is  not  unique.  A  second  asjpptotic  expansion  for 
6s  can  be  obtained  from  the  tangential  electric  field  (14)  on  s,  due 
to  a  magnetic  line  dipole  W,  and  the  impedance  boundary  condition  of 
(5a).  Carrying  out  the  necessary  calculations  yields 

GS(Q|Q')  =  6(Q|Q* )  +  Ge(Q>)  (21) 

where  the  expression  for  Ge(Q+| Q* )  can  be  obtained  from  Gh  of  (17a) 
by  simply  changing  the  a.  and_8h  to  a  and  8C.  Consequently,  incor¬ 
porating  Gs  of  (21)  intonthe  H  rield  integral  equation  of  (11a)  yields 
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identically  the  same  result  as  the  one  which  is  obtained  via  the 
E-field  Integral  equation  of  (15)  in  conjunction  with  (16b).  As 
pointed  out  earlier,  (15)  in  conjunction  with  (16b)  is  to  be  pre¬ 
ferred  over  (11a)  in  conjunction  with  (18a)  for  purposes  of  nu¬ 
merical  calculations.  It  is  interesting  to  note  that  if  one  adds 
the  zero  term,  "1-1"  to  the  integrand  of  (B46),  one  then  again 
arrives  at  (21)  after  an  appropriate  recombination  of  the  terms 
in  the  modified  integrand. 

^  It  is  observed  that  the  approximation  in  (17a)  and  (18a)  for 
Gn  (and  hence  for  Ge)  and  Gs  are  in  terms  of  tabulated  functions, 
and  they  can  therefore  be  easily  evaluated.  The  representations 
in  (17a),  (18a)  and  (16b)  are  valid  for  the  arc  lengths  Q'Q,  and 
Q'Qp  which  are  only  moderately  large.  For  larger  impedance1patch 
lengths,  i.e.,  for  t>J$t  ,  one  must  use  a  residue  series  represen¬ 
tation  of  (17b)  and  (18d)  for  Gn  and  G,  respectively.  The  de¬ 
velopment  of  the  residue  series  representations  given  in  (17b) 
and  (18b)  is  also  discussed  in  Appendix  B;  in  particular,  one  is 
referred  to  (B37)  and  (B74)  of  Appendix  B.  Such  a  residue  series 
representation  has  been  studied  extensively  by  Bremmer[21],  Logan 
[20],  and  Wai t [7] .  However,  the  alternative  representations  in 
(17a)  and  (18a)  which  are  valid  close  to  the  source  appear  to  be 
different  from  those  obtained  by  Wait [12].  Furthermore,  the  an¬ 
alysis  [9]  which  was  for  a  spherical  surface,  and  only  for  the 
TE  case,  is  not  directly  applicable  to  the  2-D  case  of  interest 
here  as  it  would  lead  to  results  that  are  not  uniformly  valid  in 
the  vicinity  of  the  source  (for  both  the  TEz  and  TMz  cases). 

The  representation  in  (17a)  and  (18a)  for  the  TE  and  TM 
cases,  respectively  consists  of  several  terms  which  have  an  im¬ 
portant  physical  significance.  The  first  set  of  terms  in  (17a) 
and  (18a)  corresponding  to  Gj;  and  G*  of  (BIO)  and  (B49)  In  Appendix 
B,  respectively  represent  the  high  frequency  currents  on  a  per¬ 
fectly  conducting  circular  cylinder;  whereas,  the  remaining  terms 
represent  a  perturbation  or  correction  arising  from  the  presence 
of  a  non-zero  surface  impedance.  Furthermore,  the  first  terms 
in  (17a)  and  (18a)  are  expressed  as  a  product  of  the  currents  on 
a  perfectly  conducting  planar  surface  and  the  perturbation  effects 
arising  from  the  finite  curvature  of  the  perfectly  conducting 
circular  cylinder.  Also,  the  second  set  of  terms  in  (17a)  and 
(18a)  which  represent  an  Impedance  perturbation  to  the  perfectly 
conducting  circular  cylinder,  have  the  same  form  as  the  impedance 
perturbation  to  a  planar  conducting  surface.  Thus,  if  these  second 
set  of  terms  were  used  in  conjunction  with  "only  those  parts"  of 
the  first  set  of  terms  which  correspond  to  the  planar  conducting 
surface,  one  would  then  obtain  the  same  asymptotic  currents  as 
those  on  a  planar  impedance  surface.  Finally,  the  last  set  of 
terms  in  (17a)  and  (18a)  are  additional  corrections  to  the  first 
two  sets  of  terms  arising  from  a  combination  of  curvature  and 
Impedance  effects.  Further  details  on  the  physical  interpretation 


of  6  for  the  completely  impedance  covered  circular  cylinder  are 
available  in  Appendix  B.  It  is  important  to  note  however,  that 
the  results  in  both  (17a)  and  (18a)  firstly  reduce  to  the  perfectly 
conducting  circular  cylinder  case  if  the  surface  impedance  van¬ 
ishes,  provided  a  slight  loss  is  introduced  in  the  external  medium 
(free  space)  before  performing  the  limit  on  (18a)  to  recover  the 
perfectly  conducting  case.  Secondly,  they  reduce  to  the  results 
for  the  planar  impedance  surface  as  the  curvature  of  the  cylinder 
becomes  vanishingly  small.  Thirdly,  the  results  in  (17a)  and  (18a) 
reduce  to  the  results  for  a  perfectly  conducting  planar  surface 
as  the  curvature  vanishes,  and  the  surface  impedance  also  vanishes. 
As  pointed  out  earlier,  the  expressions  for  the  surface  currents 
should  be  well  approximated  in  the  vicinity  of  the  source  in  order 
to  obtain  the  radiation  fields  accurately.  Since  the  radiation 
patterns  calculated  via  SC IM  are  found  to  agree  very  well  with 
those  based  on  other  trustworthy  independent  calculations,  it  is 
inferred  that  (17a)  and  (18a)  are  adequate  for  describing  the 
source  currents  in  the  vicinity  of  the  source. 

The  extension  of  the  circular  cylinder  results  in  (17a)  and 
(18a)  to  an  arbitrary  convex  cylinder  with  a  slowly  varying  curva¬ 
ture,  in  which  the  cylinder  is  completely  covered  by  an  impedance 
surface  is  carried  out  by  assuming  that  each  point  on  the  cylinder 
may  be  approximated  locally  by  a  circular  cylinder  of  the  same 
local  radius  of  curvature  at  that  point.  As  mentioned  previously, 
this  procedure  is  based  on  the  local  properties  of  high  frequency 
fields  as  employed  in  the  GTD  procedure.  Thus,  one  replaces  t 

of  (19a)  by  /9,  dt  where  Q  is  the  observation  point  of  the  currents 
on  the  impedance  patch  as  shown  in  Figure  1;  also,  the  Fock  pa¬ 
rameter,  ^  of  (19d)  is  replaced  in  the  usual  manner[22]  by  an 
integral  to  account  for  the  variation  of  the  surface  curvature. 

To  this  end,  one  proceeds  as  follows.  As  in  (19d) 

K  s  m|<t>|  =  m  |  (22a) 

where  a  is  the  radius  of  the  cylinder  and  t  is  the  arclength  defined 
in  (19c).  Since  m  (given  in  (19e))  is  assumed  to  be  slowly  varying, 
then  (22a)  can  be  differentiated,  yielding 


(22b) 
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where  p  in  (22b)  is  now  the  local  radius  of  curvature.  Integration 
of  (22b^  from  Q'  to  Q  results  with 


In  (23),  m  may  be  taken  outside  the  integral  if  one  replaces 

m  by  Jm(Q' )  m(£j)  to  represent  some  type  of  averaging  effect  to 
account  for  the  slowly  varying  curvature;  one  notices  that  such 
an  averaging  is  the  geometric  mean  of  m  between  points  Q'  and  Q. 
Furthermore,  this  averaging  effect  is  consistent  with  reciprocity. 

In  this  work,  the  length  of  the  impedance  patch  can  be  small; 
however,  it  cannot  be  too  small  such  that  the  source  at  Q1  lies 
close  to  either  edge  of  the  patch  (at  and  Q?  of  Figure  1). 

This  is  not  a  serious  restriction  sincei'n  practical  situations 
of  interest  here,  the  patch  needs  to  be  at  least  moderately  large 
enough  in  order  to  have  a  marked  effect  on  the  radiation  patterns. 
Over  such  moderate  distances  from  the  source  to  the  edges,  the 
field  component  which  continually  leaks  energy  off  the  surface 
(unlike  an  Elliot  mode  type  field  which  leaks  little  energy)  is 
in  general  sufficiently  weak  at  the  edges  of  the  patch  to  where 
the  reflection  of  this  field  which  travel  back  toward  the  source 
(while  again  shedding  energy)  may  be  neglected  in  the  calculation 
of  the  radiation  field.  Furthermore,  even  the  approximate  asymp¬ 
totic  calculation  of  these  reflections  appears  to  be  a  rather  com¬ 
plicated  task  which  is  not  considered  here.  On  the  other  hand, 
the  current  can  be  substantially  modified  by  the  reflection  of 
an  Elliot  type  mode  field  from  the  edges  of  the  impedance  patch. 
Hence,  the  present  analysis  incorporates  the  reflection  of  such 
an  Elliot  mode  field  through  the  use  of  a  local  reflection  coef¬ 
ficient  defined  at  the  edges  of  the  impedance  patch.  This  reflec¬ 
tion  coefficient  is  found  from  the  solution  to  the  canonical  problem 
of  surface  wave  reflection  from  a  discontinuity  in  planar  impedance 
boundaries  as  shown  in  Figure  5.  To  recapitulate,  one  therefore 
basically  assumes  here  that  the  total  field  incident  on  the  edges 
of  the  impedance  patch  generates  a  reflected  wave  on  the  surface 
which  is  well  approximated  by  only  the  reflected  Elliot  type  mode 
field,  since  this  mode  is  the  dominant  contributor  to  the  surface 
field.  The  remaining,  less  significant  component  of  the  reflected 
field  which  is  neglected  consists  of  propagating  and  non-propa¬ 
gating  (or  evanescent)  parts;  the  latter  part  does  not  contribute 
to  the  radiated  fields  and  it's  exclusion  is  therefore  justified. 

The  propagating  part  is  essentially  a  non-surface  wave  type  field 
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Figure  5--Geometry  of  the  canonical  problem  of  surface  wave 
diffraction  by  a  discontinuity  in  surface 
impedance  associated  with  two  part, 
planar  impedance. 


which  leaks  energy  off  the  surface;  in  general  it  contributes  suf¬ 
ficiently  weakly  to  the  radiated  field  for  impedance  patch  lengths 
which  are  not  extremely  small  and  hence  it  is  ignored. 

The  canonical  two  part  problem  in  Figure  5  above  has  been 
analyzed  via  the  Wiener-Hopf  technique  in  conjunction  with  Weinsteins 
factorization  procedure[4] .  The  latter  factorization  procedure 
leads  to  relatively  simple  Weiner-Hopf  factors.  From  [4],  it  is 
found  that  the  surface  wave  reflection  coefficient,  R  for  this 
problem  is 


n  Ji.  J  -X _  riu  « 

R  _  ^  “s  it  sinh  u  0  .  r  _  .nh-l  ®s 

Ks"JB_e0  •  *  sinh  F~ 


jrl  TMZ  case  (24) 


gTT  Jg  sinh  u 


K  =  sinh' 


;  TEZ  case 


The  plots  of  Rh  and  R.  are  illustrated  in  Figures  6  and  7,  respect¬ 
ively.  b 

h  s 

As  pointed  out  earlier,  G„’  represents  the  perturbation  to 
the  current  on  the  patch.  This  perturbation  term  arises  from  the 
truncation  of  the  impedance  surface  at  Q.  and  Q_.  In  the  following 


Figure  6— The  magnitude  and  the  phase  of  the  reflection  coefficient 

R^  for  the  TE2-case. 


Figure  7--The  magnitude  and  the  phase  of  the  reflection  coefficient 

R.  for  the  TM_  case. 
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analysis,  the  subscripts  "s"  and  "h"  are  dropped,  since  it  is  equally 
applicable  to  both  the  TEZ  and  the  TMZ  cases. 

Incorporating  the  above  values  of  R,  the  G„  can  be  calculated 
via  a  self  consistent  procedure  as  follows: 


lected  from  Q,  and 
Q2)  as  shown  7n  Fig- 
t),  can  be  found  simi- 
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The  total  incident  current,  Gt(Q.)  at  Q,  consists  of  two  terms 
which  are  the  current  G(Q.)  generatediby  the^principle  source  and 
observed  at  Q.,  and  the  cDrrent,  GB(Q,|  Q2)  re 
evaluated  at  Q,,  i.e.,  GT(Q, )  =  G(Q,)1-t-  GR(Qi 

ure  8.  The  total  incident  current  ^(Q^J^at1  Wo  '-Clll  UC  I  UUIIU  9  III! 
larly  by  interchanging  subscripts  1  4nd^2,  above.  The  term  Gfi(Q.j 
can  be  related  to  Gy(Q2)  via  the  reflection  coefficient  such  that  c 

.  -#m2  , _ 

Gr(Q,|Q2)  =  R  GT(Q?)e  ,  where  Q7Q7  denotes  the  arclength 

between^!),  and  (L.  Similarly,  the  tefnrGR(QJ  Q.)  (which  arises 
in  the  process  ot  calculating  G,(Q?))  can'be^written  in  terms  of 
R,  GT(Q.)  and  (L0,  •  One  may  th4n  Solve  for  GT(Q.)  and  GT(Q9). 

It  can  be  shown^tnat:  _  1 

G(Q,)+R  e  W  G(Q?) 

i - .  (26) 
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Figure  8— Incident  and  reflected  rays  associated  with 
the  self  consistent  method. 


Interchanging  subscripts  1  and  2  in  (26)  gives  the  corresponding 
expression  for  Gt(q^).  gne  notes  that  (26)  reduces  to  a  more  familiar 
result  when  the  source  point  is  located  at  the  center  of  the  patch 


■JPQiQ?  i 

so  that  6(0,)  =  G(Q2);  in  which  case  G, (Q, )  =  G(Q,)  (1-R  e  )  . 

The  perturbation  current  at  a  point  Q  on  the  paten  is  then  given 
by  the  following  equation 


-JBQi  Q 

Gr(Q)  =  R  Gt(Qj)  e  1 


+  RGt(Q2)  e 


(27) 


As  the  value  of  Z  /Z  becomes  smaller  than  0.5,  the  correspond¬ 
ing  values  of  |R  |  becomes0 larger  than  0.9  as  seen  from  Figure 
7.  Such  a  largesvalue  of  the  reflection  coefficient  R  implies 
via  (27)  that  a  very  large  number  of  multiple  bounces  Between  the 
edges  of  the  impedance  patch  remain  significant  in  the  calculation 
of  the  TM  Elliot  mode  field.  In  such  a  situation,  the  accuracy 
of  the  present  solution  (SCIM)  may  be  impaired,  because,  a  small 
error  in  the  approximation  to  the  Elliot  mode  field  which  Is  di¬ 
rectly  Incident  from  the  source  onto  the  edges  of  the  patch  could 
be  compounded  through  each  subsequent  reflection  (bounce)  from 
the  edges.  In  practice,  one  would  not  employ  values  of  |Z  /Z  | 
for  which  |R  |  is  large  anyway,  as  this  is  obviously  not  desirable 
from  an  antenna  design  point  of  view. 

Furthermore,  the  earlier  assumption;  that  is,  the  propagating, 
non-surface  wave  type  fields  which  are  generated  by  the  "reflections" 
of  Elliot  mode  from  the  edges  contribute  sufficiently  weakly  to 
the  radiated  field,  does  not  hold  for  such  large  values  of  the 
reflection  coefficient.  Excluding  this  non-surface  wave  field  results 
in  an  inaccurate  representation  of  the  currents  on  the  patch  in 
the  vicinity  of  the  edges.  However,  as  indicated  earlier,  an  ac¬ 
curate  representation  of  the  currents  on  the  impedance  patch  is 
essential  for  accuracy  in  the  SCIM.  One  notes  that  in  the  SCIM, 
an  Integration  of  the  component  of  the  current  reflected  from  the 
edges  gives  rise  to  fields  diffracted  not  only  from  the  edges  on 
which  the  reflected  component  terminates,  but  also  from  the  op¬ 
posite  edges  where  the  reflected  component  originates.  The  latter 
"spurious  diffracted  field"  contribution  is  sensitive  to  errors 
in  the  currents  at  the  edges;  however,  these  errors  are  generally 
not  large  if  the  reflection  coefficient  is  sufficiently  small. 

An  alternative  approach  for  Including  the  effects  of  multiple  re¬ 
flections  in  the  total  radiation  field  calculation  is  to  simply 
superpose  the  solution  based  on  the  SCIM  but  without  including 
the  correction  currents  resulting  from  the  reflections  of  the  Elliot 
modes,  with  the  part  of  the  GTD  ray  solution  of  [4]  for  calculating 
the  diffracted  fields  which  result  only  from  these  reflections. 

The  latter  GTD  diffracted  field  contribution  is  calculated  via 
the  self  consistent  approach  for  summing  up  the  multiple  reflec¬ 
tions  of  the  Elliot  modes  as  indicated  earlier.  It  is  noted  that. 
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the  GTD  solution  does  not  require  an  exact  knowledge  of  the  form 
of  the  currents  In  the  vicinity  of  the  edges,  hence  it  is  suitable 
for  the  evaluation  of  the  contribution  to  the  radiation  pattern 
from  the  multiple  reflections.  Furthermore,  this  alternative  pro¬ 
cedure  is  attractive  from  the  point  of  making  the  numerical  radi¬ 
ation  field  calculations  more  efficient. 

C.  Numerical  Evaluation  of  the  Surface  Currents 

on  the  Impedance  Patch 

The  last  step  of  this  analysis  is  to  evaluate  the  integrals 
in  (11a),  (lib)  and  (15),  numerically  over  the  impedance  patch. 

As  pointed  out  earlier,  H® *m( Ql P * )  is  the  maqnetic  field  on  the 
impedance  patch  due  to  an°electric  (magnetic)  line  source  at  P'; 

Hq  ( Q | P 1 )  is  a  well  behaved  function  of  the  arclength  on  the 

cylinder  and  suitable  for  a  numerical  integration  technique.  The 
explicit  form  of  this  magnetic  field  is  given  in  [1,2]  and  will 
not  be  repeated  here. 

The  term  nxTT.  ( Q | Q ' )  in  (lib)  is  related  to  (17a)  via  (16a). 

As  it  can  be  seenDfrom  (17a),  this  term  has  a  log-type  singularity 
(small  argument  approximation  for  H^z'(:<t))  as  the  variable  of 
Integration  (i.e.,  the  arclength)  goes  through  the  source.  Such 
integrable  singularities  can  be  treated  numerically  in  the  computer 
via  Gaus'  method.  Thus,  the  TE  case  presents  no  numerical  dif¬ 
ficulties.  However,  the  term  nxH.  (Q|Q‘)  in  (11a)  is  much  too 
singular  for  a  numerical  integration  scheme.  In  fact,  from  (18a), 
with  the  useage  of  the  small  argument  approximation  of  Hjz(kt), 
one  can  show  that  this  singularity  is  in  the  fornuof  an  inverse 
square  of  the  arclength  type  (i.e.,  l/( arc  length)^).  Such  singu¬ 
larities  are  not  suitable  for  a  numerical  integration  procedure 
and  therefore  one  must  reduce  the  singularity  by  means  of  the 
distributed  source  technique  before  any  calculation  work  can  be 
done.  To  this  end,  one  proceeds  as  follows  for  the  nxTT.  of  the 
TM  case.  Let  the  source  be  a  distributed  sinusoid  sucn  that  the 
distribution  function,  A(x')  is  given  by 


Afv'l  -  sin  k(w-lx'l) 
A{x  '  sin  kw 


(28) 


where  2w  is  the  total  extent  of  the  source.  The  resulting  current 
on  the  patch  from  the  slnsoidally  distributed  source  can  be  found 
by  weighting  the  currents  generated  by  an  infinitesimal  source 
with  the  distribution  function,  A(x’)  and  then  Integrating  it  from 
-w  to  w,  i.e., 
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(29) 


nxHb(Q)  = 


/  A(x')  nxHb(Q(Q')  dx'  , 


-w 


where  Hb(Q)  and  H.(Q|Q')  In  (29)  are  the  magnetic  fields  observed 
at  a  point  0  on  tne  surface  of  the  patch  due  to  a  sinusoidally 
distributed  magnetic  line  dipole,  and  an  infinitesimal  magnetic 
line  dipole,  respectively. 

As  mentioned  above,  the  singularity  of  nxHb(Q | Q ’ )  in  (29) 
has  the  Hj^'(kt)/kt  form  (see  Equation  (18a))  wnich  can  be  rewritten 
as 


Hx(kt) 

ict~ 


H^2)(kt) 


(30) 


In  general,  if  a  function,  F( k I x-x * ( )  can  be  rewritten  in 
terms  of  another  function  f(kjx-x'|)  via 


F(k|x-x'|)  -  fl  +4 - yl  f(k|x-x'|)  (31) 

d2(k|x-x‘|)2  J 


then  the  integral  of  F (k ( x-x ' I )  with  A(x')  over  the  interval  (-w,w) 
can  be  expressed  in  terms  of  f(k|x-x'|)  as  follows 


I  = 


w 

/ 


-w 


A(x')  F(k jx-x' | ) dx  * 


/  A(x' )  fl  +  — - 7V(k|x-x'  ( ) dx ' . 

-w  V  d(k|x-x'|)v 

(32) 


Repeated  application  of  integration  by  parts  to  (32)  yields 


I  =  A(x')  f ’ (k j x-x ' | )  |  -  A'(x')  f(kjx-x'j)  |  + 


-w 


/  dx' |A(x' ) 
-w 


+A”(x ' ) | f (k ) x-x 1 | )  ,  (33) 
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where  primes  over  the  functions  indicate  derivatives  with  respect  to 
the  argument. 

The  first  term  in  (33)  is  zero  because  A(±w)=0.  Furthermore  the 
last  term  in  (33)  is  also  zero,  because  A"(x')  =-A(x').  Hence  (33) 
becomes 

I  -  -A'(x')f(k|x-x'|)|°  -  A'(x')f(k|x-x'|)|J  .  (34) 

-w 

From  the  definition  of  A(x'),  it  can  be  shown  that 


3F  A(x,) 


+  cos 


(35) 


Substituting  (35)  into  (34)  yields  the  desired  relationship 


b  /  j2 

/  A(x')  (  1  +4 - 

-b  \  d(k|x 


— — 2^f(k|x-x*|)dx‘  s-^r^Cf(klx"wl) 


+  f(k|x+w|)-2  cos  kw  f(k|x|)] 


Substituting  (30)  into  (18a)  and  employing  (36)  yields 


Gs 


H02)(kt1>Fl(5,>  -  Ij  H<2)(kt,){l-2zieZl 


FltklVflidt’  +  fsA(t0)c, 


(37) 


where  5s  is  the  integral  of  Gs  of  (18a)  weighted  by  A(t')  over 
the  interval  -w  to  w.  The  t  ,,  t  and  t<,  are  the  arc lengths  be¬ 
tween  the  observation  point  Iron  the  patch  and  the  points  -w,  o, 
and  w  of  the  distributed  source,  respectively  as  shown  in  Figure 


9.  Sj  and  z^  are  defined  in  terms  of  tj  in  a  similar  fashion, 
t  <|  „ 

such  that  £.j  =  m  —  and  z^  =  jbg  kt^.  The  result  in  (37)  can 


be  generalized  to  an  arbitrary  convex  cylinder  with  a  slowly  vary¬ 
ing  curvature  via  (23). 
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Figure  9— Sinusoidally  distributed  line  source 
and  related  parameters. 


The  resulting  current  on  the  patch  due  to  a  distributed  source* 
A(t')  can  now  be  written  as  follows 


nxHb(Q|Q-)  =  -jYflMd(GS  +  G|)z 


(38) 


~  c  e 

Gp  can  be  obtained  from  (27),  if  one  uses  G  instead  of  G  . 

As  it  can  be  seen  from  (37),  the  highest  singularity  in  (38) 
is  a  log- type  singularity;  hence,  the  resulting  integrals  upon 
substitution  of  (38)  into  (lib)  can  be  evaluated  numerically  as 
principal  values. 

The  alternative  formulation  of  the  TM  discussed  above  In 
Section  I IB  is  recommended  for  the  following  reasons.  Firstly,  the 
currents  of  (18a)  to  be  employed  in  the  reaction  integral  formula¬ 
tion  given  in  (11a)  are  highly  singular  at  the  source  and  therefore 
it  is  very  difficult  to  obtain  a  satisfactory  mathematical  repre¬ 
sentation  for  the  fields.  Secondly  an  attempt  at  smoothing  out 
these  singularities  (as  it  is  done  above)  for  the  purposes  of 
numerical  treatment  leads  to  a  rather  complicated  and  slow  com¬ 
puter  code.  It  is  however  important  to  point  out  that  the  results 
obtained  from  (11a)  agree  very  well  with  other  independent  checks 
such  as  MM-GTD  and  GTD  techniques.  Nevertheless,  the  electric 
field  integral  representation  given  in  (15)  for  the  TM  case  is 
superior  to  the  one  given  in  (11a)  because  the  formulation  in  (15) 
leads  to  source  singularities  which  can  be  numerically  Integrated 
in  the  principle  value  sense;  this  numerical  treatment  is  not  com¬ 
plicated,  and  It  leads  to  faster  computer  codes.  Furthermore, 
the  computer  code  for  numerically  processing  the  electric  field 
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representation  of  (15)  for  the  TM  case  can  be  obtained  from  the 
code  for  the  TE  case  (lib)  with  only  minor  changes;  hence,  both 
cases  can  be  incorporated  into  a  single  computer  code.  It  is  how¬ 
ever  important  to  note  that.  In  the  process  of  numerical  integration, 
one  should  divide  the  impedance  patch  into  small  enough  segments 
so  that  the  integral  subroutines  can  sample  the  variations  of  the 
currents  on  the  patch  accurately.  It  is  observed  that  the  segments 
which  are  as  long  as  half  an  electric  wavelength;  i.e.,  w/B  can 
accurately  be  Integrated  via  a  5-point  Gaussian  quadrature  formula. 
The  B  (8h  for  TE  and  8,  for  TM  )  is  given  in  (19g)  and  (20d). 


CHAPTER  III 

RELIABILITY  OF  THE  SURFACE  CURRENT  INTEGRATION  METHOD  (SCIM) 


This  chapter  is  designed  to  establish  the  validity  of  the 
present  analysis  by  comparing  it  with  other  independent  methods. 
These  numerical  comparisons  of  the  SCIM  are  illustrated  only  for 
the  circular  cylinder  geometry  of  Figure  10.  Nevertheless,  as 
pointed  out  earlier,  the  present  analysis  based  on  the  SCIM  can 
also  be  directly  employed  to  treat  convex  cylinders  of  arbitrary 
cross  sections.  For  convenience,  the  source  is  located  in  the 
middle  of  the  impedance  patch.  The  impedance  patch  width  is  de¬ 
fined  in  terms  of  the  angle,  measured  from  the  source.  The 
angles,  <j>  and  4>n  specify  the  points  on  the  impedance  patch  and 
the  observation*direction  in  the  far  field,  respectively. 

The  SCIM  is  compared  with  three  independent  methods.  Two 
of  these  methods  are  based  on  the  different  mathematical  models 
of  the  problem,  and  they  are  called  the  MM-GTD  and  the  GTD  formu¬ 
lations.  Some  experimental  results  are  also  included  as  a  third 
method.  It  is  noted  that  like  the  SCIM,  the  MM-GTD  and  GTD  formu¬ 
lations  are  also  asymptotic  in  the  sense  that  they  are  developed 
for  electrically  large  cylinders.  The  MM-GTD  and  the  GTD  formula¬ 
tions  are  discussed  in  detail  elsewhere  [2,4]  and  some  of  the  main 
features  of  these  methods  have  been  discussed  in  Chapter  I,  hence 
only  a  short  summary  of  these  methods  will  be  repeated  here  for 
completeness. 

The  MM-GTD,  a  hybrid  technique,  combines  the  method  of  moments 
and  the  geometrical  theory  of  diffraction.  An  application  of  the 
generalized  reaction  theorem,  in  a  manner  analogous  to  that  employed 
in  the  development  of  the  result  in  Figure  11,  yields  the  governing 
integral  equation  employed  in  the  present  MM-GTD  formulation  as: 


/  VQ|Q,)  *  fy$)nxnxHj(q|P)]dA  =-flJJ(Q'|P)  -z+Hoz(P|Q')  .  z 
$1 

For  TEZ  case  (39a) 


•  fZs(Q)nxnxTf^(Q|P)]<Ui  =-T^(Q'|p)  •  f-r^PlQ')  • 

I 

For  TMz  case.  (39b) 
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In  (39),  fl  and  I  are  produced  by  the  source  WT=z6( |p-pn,  | ) ;  whereas 
TT  and  IT  02 are  produced  by  M  =t6( |p-pn, | ) .  As  in  the  SClH,  the  sub¬ 
scripts  "8*  and  "b"  in  (39),  denotes  the  fields  of  the  unperturbed  and 
perturbed  problem  respectively.  Equations  (39a)  and  (39b),  are  first 
order  inhomogeneous  integral  equations  with  unknown  "FT  and  TTr  re¬ 
spectively,  and  they  can  be  solved  via  the  moment  method  technique. 

What  differentiates  this  method  from  the  conventional  moment  method 
solution  of  integral  equations  pertaining  to  electromagnetic  prob¬ 
lems  Is  the  fact  that  in  the  former  one,  the  free  space  Green's 
function  is  used  as  the  kernel  of  the  integral  equation;  whereas 
in  the  present  analysis,  a  uniform  GTD  approximation  for  the  per¬ 
fectly  conducting  convex  cylinder  Green's  function  constitutes 
the  kernel.  This  is  the  reason  behind  choosing  the  name  MM-GTD. 

At  the  present  time,  the  MM-GTD  solution  for  the  TE  case  appears 
to  be  very  accurate;  on  the  other  hand,  the  TM  casi  is  expected 
to  be  only  reasonably  accurate  since  further  work  is  needed  to 
refine  the  MM-GTD  procedure  for  the  TMz  case  as  indicated  in  [5]. 

The  GTD  method  is  a  ray  optical  method,  and  the  ray  picture 
for  this  problem  is  illustrated  in  Figure  2.  The  direct  (or  geo¬ 
metrical  optics)  ray  field  and  the  fields  diffracted  from  the  edges 
of  the  patch  are  calculated  from  appropriate  canonical  problems 
for  TM  and  TE  cases  as  explained  in  [4].  It  is  noted  that  the 
GTD  solution  presented  in  [4]  is  restricted  to  only  those  impedance 
surfaces  which  can  support  an  Elliot  type  surface  ray  mode  field. 

As  indicated  above,  an  experimental  verification  of  the  ac¬ 
curacy  of  the  SCIM  solution  is  also  carried  out  in  this  study. 

This  verification  is  for  the  TE  case.  The  experimental  set-up 
consists  of  a  5  ft  long,  3*s  incn  radius  cylinder,  which  models 
the  infinitely  long  perfectly  conducting  cylinder.  An  x-band  wave¬ 
guide  feeds  an  axial  slot  in  this  cylinder  as  illustrated  in  Figure 
11.  Polyethylene  layers  which  are  4  inches  wide,  and  0.058  inches 
thick  are  employed  to  simulate  the  impedance  patch.  The  value 
of  the  impedance  can  be  changed  by  altering  the  number  of  poly- 
thelene  layers.  The  measurements  are  conducted  at  10.77  GHz. 

At  this  operating  frequency,  the  complex,  dielectric  constant  of 
polyethylene  is  found  to  be  e  =  (2.25  +  j9xl0je  ;  e  being  the 
permittivity  of  the  free  space.  0  0 

Comparisons  of  the  results  for  the  radiation  patterns  based 
on  the  SCIM  with  those  based  on  the  MM-GTD,  the  GTD,  and  the  experi¬ 
mental  data  are  presented  in  Figures  12  through  28. 

The  GTD  patterns  presented  here  appear  to  be  reasonably  accurate. 
However,  it  is  observed  that  the  GTD  results  in  [4]  for  the  circular 
cylinder  TMZ  case  in  general  need  to  be  refined  further  to  make  them 
as  reliable  as  for  the  TE  case,  especially  for  a  wider  range  of 
parameters  than  those  considered  in  the  GTD  calculations  presented 
here.  In  the  case  of  a  cylinder  with  a  very  small  curvature 


Figure  12— Radiation  patterns  of  a  magnetic  line  source  on 
perfectly-conducting  planar  and  circular  cylindrical 
surfaces  which  are  partly  covered  with  an 
impedance  surface  patch. 


(ka=360),  the  asymptotic  analysis  (SCIM)  in  Figures  12,  13,  14, 

25  and  26  are  also  compared  with  the  results  obtained  from  the  GTD 
analyses  for  an  impedance  patch  on  a  perfectly  conducting  ground 
plane  (which  is  the  limiting  case  of  a  cylinder  with  vanishing 
curvature.)  These  radiation  patterns  of  an  antenna  on  a  flat 
surface  can  be  used  to  predict  with  reasonable  accuracy  the  perti¬ 
nent  parameters  of  the  radiation  patterns  in  the  lit  region  of 
a  similar  antenna  located  on  a  cylindrical  surface  with  small  curva¬ 
ture. 


In  Figures  15  through  22,  and  27  and  28,  the  SCIM  is  compared 
with  the  GTD  and  MM-GTD  results.  It  is  observed  that  the  results 
obtained  from  these  three  independent  methods  compare  very  well 
with  each  other.  However,  it  should  be  pointed  out  that  there 
are  no  MM-GTD  results  presented  for  the  TM  case.  This  is  because 
it  is  very  cumbersome  to  obtain  a  satisfactory  MM-GTD  solution 
for  this  case  at  the  present  time,  as  has  been  mentioned  before. 
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Figure  13.  Radiation  patterns  of  a  magnetic  line  source  on 
perfectly-conducting  planar  and  circular  cylindrical 
surfaces  which  are  partly  covered  with  an 
impedance  surface  patch. 


The  relative  merits  of  these  three  methods  (i.e.,  SCIM,  MM- 
GTD  and  GTD)  are  presented  in  Table  I.  In  Figures  22 ,  23,  and 
24,  where  the  SCIM  is  compared  with  experimental  results,  it  is 
modified  by  taking  the  impedance  surface  as  being  slightly  lossy 
to  account  for  the  loss  in  polyethelene  panels.  The  loss  is  in¬ 
troduced  into  ..he  solution  as  an  exponential  decay  along  the  im¬ 
pedance  surface.  From  this  point  of  view  the  solution  may  be  con¬ 
sidered  as  a  perturbation  to  the  lossless  case.  The  SCIM  is  further 
modified  by  introducing  the  effect  of  the  finite  source  distribution 
in  the  axial  slot  via  a  multiplicative  slot  pattern  factor.  The 
polyethylene  layers  are  modeled  as  an  impedance  surface.  These 
polyethylene  layers  which  are  0.058  inches  and  0.106  inches  thick 
for  the  two  experimental  cases  may  be  approximated  by  an  equiva¬ 
lent  impedance  of  Z  =j70fi  and  Z  =jl50fl,  respectively;  this  equiva¬ 
lent  inductive  impedance  is  based  on  the  flat  surface  approximation 
and  the  corrections  to  this  impedance  arising  from  the  curvature 
of  the  cylinder  have  been  ignored  since  the  cylinder  is  over  twenty 
wavelengths  in  circumference. 
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Figure  14.  Radiation  patterns  of  a  magnetic  line  source  on 
perfectly-conducting  planar  and  circular  cylindrical 
surfaces  which  are  partly  covered  with  an 
impedance  surface  patch. 


Finally,  it  is  recognized  that  the  TM  case  is  very  nearly 
a  dual  of  the  TE  configuration.  The  principle  difference  is  that 
the  TM  source  is  a  magnetic  line  dipole  source  which  is  of  course 
not  a  dual  of  the  magnetic  line  source  and  the  conducting  part 
of  the  surface  in  the  TE  case  is  not  replaced  by  it’s  dual  for 
the  TM  case.  Nevertheless,  by  a  proper  choice  of  surface  im¬ 
pedances,  it  is  possible  to  generate  patterns  (for  the  TM  case) 
similar  in  nature  to  the  smooth  patterns  of  the  TE  case  presented 
in  this  chapter.  The  following  results  for  the  TMZ  case  are  taken 
for  the  range  of  impedances  such  that  high  rippleszoccur  in  the 
radiation  pattern.  This  is  done  merely  to  compare  the  SCIM  with 
the  GTD  solution  in  the  domain  where  they  both  are  valid,  and  to 
supplement  the  computations  for  the  smooth  patterns  for  the  TE 
case.  It  is  emphasized  that  smooth  patterns  can  also  be  generated 
for  the  TM  case  as  is  shown  in  Figures  46  through  51. 
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Figure  15--Radiation  patterns  of  a  magnetic  line  source  on 
perfectly-conducting  circular  cylinder  which  is 
partly  covered  with  an  impedance  surface  patch. 
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Figure  16--Radiation  patterns  of  a  magnetic  line  source  on 
perfectly-conducting  circular  cylinder  which  is 
partly  covered  with  an  impedance  surface  patch. 
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Figure  17--Radiation  patterns  of  a  magnetic  line  source  on 
perfectly-conducting  circular  cylinder  which  is 
partly  covered  with  an  impedance  surface  patch. 
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Figure  18— Radiation  patterns  of  a  magnetic  line  source  on 
perfectly-conducting  circular  cylinder  which  is 
partly  covered  with  an  impedance  surface  patch. 
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Figure  19--Radiation  patterns  of  a  magnetic  line  source  on 
perfectly-conducting  circular  cylinder  which  is 
partly  covered  with  an  impedance  surface  patch. 
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Figure  20— Radiation  patterns  of  a  magnetic  line  source  on 
perfectly-conducting  circular  cylinder  which  is 
partly  covered  with  an  impedance  surface  patch. 
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Figure  21--Radiation  patterns  of  a  magnetic  line  source  on 
perfectly-conducting  circular  cylinder  which  is 
partly  covered  with  an  impedance  surface  patch. 
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Figure  22— Radiation  patterns  of  a  magnetic  line  source  on 
a  perfectly-conducting  cylinder. 
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Figure  23— Radiation  patterns  of  a  magnetic  line  source  on 
perfectly-conducting  circular  cylinder  which  is 
partly  covered  with  an  impedance  surface  patch. 
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Figure  24--Radiation  patterns  of  a  magnetic  line  source  on 
perfectly-conducting  circular  cylinder  which  is 
partly  covered  with  an  impedance  surface  patch. 
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Figure  25--Radiation  patterns  of  a  magnetic  line  dipole  source 
on  perfectly-conducting  planar  and  circular  cylindrical 
surfaces  which  are  partly  covered  with  an 
impedance  surface  patch. 
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Figure  26--Radiation  patterns  of  a  magnetic  line  dipole  source 
on  perfectly-conducting  planar  and  circular  cylindrical 
surfaces  which  are  partly  covered  with  an 
impedance  surface  patch. 
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Figure  27--Radiation  patterns  of  a  magnetic  line  dipole  source 
on  perfectly-conducting  circular  cylinder  which  is  partly 
covered  with  an  impedance  surface  patch. 


Figure  28--Radiation  patterns  of  a  magnetic  line  dipole  source 
on  perfectly-conducting  circular  cylinder  which  is  partly 
covered  with  an  impedance  surface  patch. 
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It  is  observed  from  these  comparisons  that  the  asymptotic 
solution  (SCIM)  indeed  agrees  very  well  with  the  other  independent 
methods  of  analysis. 

In  the  following  paragraphs,  the  bounds  for  the  asymptotic 
analysis  (SCIM)  are  studied.  The  basic  parameters  of  interest 
occuring  in  the  SCIM  are:  (a)  the  distance  from  the  source  to 
the  edge  of  the  patch,  (b)  the  radius  of  the  cylinder,  and  (c) 
the  value  of  the  surface  impedance.  The  range  of  these  parameters 
over  which  the  SCIM  remains  accurate  have  been  established  in  an 
approximate  fashion  after  an  extensive  numerical  study.  In  par¬ 
ticular,  for  a  given  value  of  impedance  and  the  cylinder  radius, 
there  is  a  maximum  patch  size  beyond  which  the  asymptotic  repre¬ 
sentation  used  in  the  SCIM  for  the  currents  on  the  impedance  sur¬ 
face  are  strictly  not  valid.  As  indicated  previously  in  Chapter 
II,  the  surface  current  on  the  impedance  patch  is  assumed  to  consist 
partly  of  the  Elliot  mode  or  surface  wave  type  component  which 
propagates  from  the  source  to  the  edge  of  the  patch,  and  secondly 
it  consists  of  the  component  which  is  termed  as  the  non-surface 
wave  component,  since  the  latter  is  the  only  current  which  propa¬ 
gates  to  the  edge  of  the  patch  if  the  impedance  surface  can  not 
support  an  Elliot  mode  field.  Thirdly,  there  is  a  reflected  com¬ 
ponent  of  the  current  to  account  for  the  reflections  of  the  current 
incident  on  the  edges  of  the  patch.  Also,  as  pointed  out  earlier 
this  reflected  component  of  the  current  is  assumed  to  be  produced 
by  the  Elliot  mode  field  or  surface  wave  component  which  inpinges 
on  the  edge  of  the  patch.  A  few  typical  plots  of  the  magnitudes 
of  the  equivalent  magnetic  surface  current  densities  which  exist 
on  an  impedance  patch  located  on  a  perfeclty  conducting  circular 
cylinder  are  illustrated  in  Figures  29  through  34.  The  ^urr^nt 
densities  (Z  nxnxH ( Q | Q 1 )  of  (16a)  for  the  TE  -case  and  nxE.  (Q  ) Q ’ ) 
of  (16b)  forthe  TM  -case)  are  plotted  as  a  function  of  the  arc- 
length,  t=a 1 4> | /X ,  measured  from  the  source  placed  at  the  center 
of  the  impedance  patch  where  X=2-rr/k  and  |<j>|  is  the  angle  between 
the  source  point  and  the  observation  point  on  the  impedance  patch. 

As  can  be  seen  from  these  typical  current  plots,  the  amplitude 
of  the  non-surface  wave  component  of  the  current  begins  to  grow 
beyond  a  certain  maximum  distance,  t  /2  from  the  source.  The  non¬ 
surface  wave  current  approximation  i?  not  expected  to  be  accurate 
beyond  the  immediate  neighborhood  of  this  distance,  t  /2.  For 
the  TE  -case,  an  empirical  equation  can  be  found  relaxing  the  maxi¬ 
mum  total  patch  size,  t  (as  defined  above  and  measured  in  wave¬ 
lengths;  i.e.,  2a|<j>  |/A;  to  the  magnitude  of  the  normalized  im¬ 
pedance,  |z  I  =  |zJf/Z  anc*  tlie  normalized  radius  of  the  cylinder, 
ka.  For  the  asymptotic  forms  of  the  currents  which  are  given  in 
section  I  IB  to  order  l/(ka)  ,  such  an  emperical  equation  has  the 
following  form: 

ka  =  [|zs|  +  0.04]  (3.1  t*  -  2.6  t  +  2.8)  +  3.6.  (40) 
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Figure  29— Magnitude  of  the  equivalent  magnetic  current  density 
on  an  impedance  patch  {which  partly  covers  a  circular 
cylinder)  as  a  function  of  the  arclength  measured 
from  the  source.  R. =0.071/14°. 
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Figure  30— Magnitude  of  the  equivalent  magnetic  current  density 
on  an  impedance  patch  (which  partly  covers  a  circular 
cylinder)  as  a  function  of  the  arclength  measured 
from  the  source.  R.  =0.2173/24°. 
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Figure  31--Magnitude  of  the  equivalent  magnetic  current  density 
on  an  impedance  patch  (which  partly  covers  a  circular 
cylinder)  as  a  function  of  the  arclength  measured 
from  the  source.  Rh=0. 2173/24°. 


Figure  32--Magnitude  of  the  equivalent  magnetic  current  density 
on  an  impedance  patch  (which  partly  covers  a  circular 
cylinder)  as  a  function  of  the  arclength  measured 
from  the  source.  R  =0.7825/121°. 
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Figure  33--Magnitude  of  the  equivalent  magnetic  current  density 
on  an  impedance  patch  (which  partly  covers  a  circular 
cylinder)  as  a  function  of  the  arclength  measured 
from  the  source.  R  =0.7825/121°. 
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Figure  34--Magnitude  of  the  equivalent  magnetic  current  density 
on  an  impedance  patch  (which  partly  covers  a  circular 
cylinder)  as  a  function  of  the  arclength  measured 
from  the  source.  R  =0.9292/129°. 
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It  is  noted  that,  (40)  is  a  very  conservative  upper  bound 
for  both  the  TE  and  the  TM  cases.  This  is  particularly  true 
for  the  type  ofzimpedances  Which  can  allow  an  Elliot  type  of  surface 
wave  mode  to  exist,  since  this  is  the  dominant  component  of  the 
current  and  furthermore  it  remains  valid  even  for  distances  larger 
than  t  /2.  Figures,  presented  in  (12)  through  (28),  indicate  that 
it  is  possible  to  obtain  good  agreement  between  the  SCIM  and  other 
independent  methods  for  patch  sizes  well  beyond  those  allowed  by 
(40).  Typically,  it  is  in  general  possible  to  treat  half  patch 
widths  of  the  order  of  t  with  good  accuracy. 

For  impedance  patches  which  are  larger  than  the  limiting  size 
indicated  above,  one  should  use  the  asymptotic  form  of  the  currents 
given  in  (17a)  and  (18a)  up  to  t  and  then  switch  to  the  creeping 
wave  representation  for  the  currents  which  is  explained  in  parts 
(IC)  and  (IIC)  of  Appendix  B,  for  the  TE  and  TM  cases,  respec¬ 
tively.  2  2 


The  component  of  the  current  reflected  from  the  edges  of  the 
impedance  patch  is  calculated  by  simply  multiplying  the  incident 
Elliot  mode  field  with  a  surface  wave  reflection  coefficient  R. 


which  characterizes  the  discontinuity  in  the  impedance  at  the  edges. 
This  quantity,  R.  is  found  in  [4]  via  the  solution  to  the  problem 


of  the  surface  wave  diffraction  by  a  two  part  planar,  impedance 
surface.  Multiple  reflections  between  the  edges  of  the  patch  are 
then  summed  in  a  self  consistent  manner  as  explained  in  section 
I IB .  Plots  of  R.  (corresponding  to  TE  case)  and  R  (corresponding 
to  TM  case)  are  illustrated  in  Figures  6  and  7,  respectively  as 
a  function  of  the  surface  reactance  which  is  normalized  to  the 
free  space  impedance.  It  is  noted  from  these  plots  that  R  is 
significantly  higher  than  R,  for  values  of  the  normalized  react¬ 
ance  which  are  less  than  unity.  Consequently,  the  component  of 
the  current  reflected  from  the  edge*  of  the  impedance  patch  is 
expected  to  play  a  significant  role  in  the  radiation  pattern  calcu¬ 
lations  for  the  TM  case. 


In  general,  the  SCIM  is  not  expected  to  be  accurate  for  the 
case  when  the  magnitude  of  the  reflection  coefficient  |R  |  approaches 

B 

unity.  However,  as  mentioned  at  the  end  of  section  I  IB ,  this  is 
not  a  serious  restriction  on  the  SCIM  since  an  antenna  wi^h  such 
a  large  reflection  coefficient  is  not  desirable.  Linder  the  present 
assumption,  the  dominant  reflected  component  of  the  current  is 
produced  by  the  incident  Elliot  mode  field.  In  order  for  this 
assumption  to  hold,  one  requires  that  the  size  of  the  impedance 
patch  be  sufficiently  large.  From  an  extensive  numerical  study, 
it  has  been  approximately  determined  that  this  minimum  impedance 
patch  length  for  the  case  when  the  impedance  surface  supports  an 
Elliot  mode  field  is  n/2B  where  3  for  both  the  TE  and 
TM  cases  defined  in  (19g)  and  (20d)  respectively.  This  study 


basically  indicates  how  far  the  edge  of  the  patch  has  to  be  to 
lie  outside  the  "effective  source  regions"  where  the  component 
of  the  currents  propagating  away  from  the  source  are  usually  large 
(due  to  the  singularity  at  the  source).  It  is  noted  that  tt/28 
is  generally  very  small  compared  to  t  so  that  the  current  approxi 
mation  employed  in  SCIM  is  indeed  valTd. 


t 


CHAPTER  IV 

RELATION  OF  RADIATION  PATTERN  PROPERTIES  TO  THE  PERTINENT 
PARAMETERS  OF  THE  PROBLEM 


In  this  chapter  the  dependence  of  the  radiation  patterns  on 
certain  parameters  is  discussed.  It  is  observed  that  the  number 
of  lobes  in  the  radiation  pattern  within  the  lit  region  is  a  function 
of  the  width  of  the  impedance  patch  and  the  value  of  the  impedance. 
For  the  TM  case  with  a  capacitive  impedance  patch,  this  number 
is  primarily  governed  by  the  width  of  the  impedance  patch  measured 
in  wavelength,  x  (where  as  before  X=2rr/k).  A  simple  ray  picture 
of  the  GTD  may  be  employed  to  explain  why  the  number  of  the  lobes 
in  the  lit  region  depends  primarily  on  the  impedance  patch  length. 

In  this  picture,  the  principle  source  and  the  edges  of  the  impedance 
patch  may  be  considered  as  three  line  sources.  For  the  TM  case, 
these  line  sources  have  similar  pattern  factors.  Hence  a  simple 
array  analysis  indicates  that  the  number  of  lobes  in  the  radiation 
pattern  within  the  lit  region  is  approximately  given  by  the  fol¬ 
lowing  table. 


Table  2 


Width  of  the  impedance 
patch  (in  X) 

Approximate  number  of 
lobes  in  the  TM  radi¬ 
ation  pattern  (Tit  region) 

Less  than  0.5 

1 

0.5-1. 5 

3 

1.5-2. 5 

5 

2. 5-3. 5 

7 

3. 5-4. 5 

9 

Table  2  is  in  full  agreement  with  the  observed  radiation  pat¬ 
terns  for  the  TM  -case.  Having  a  peak  in  the  broadside  direction 
(i.e.,  4>n=0°)  isza  characteristic  feature  of  this  case. 

For  the  TE  case,  with  an  inductive  impedance  patch,  it  is 

difficult  to  relate  the  number  of  lobes  in  the  lit  region  to  just 

the  impedance  patch  length  via  a  simple  array  analysis  because 
the  pattern  factor  of  the  principle  source  located  at  the  center 

of  the  patch  is  different  from  the  pattern  factors  of  the  equiva¬ 

lent  sources  of  the  fields  diffracted  from  the  edges  at  and 
Q-.  In  general  for  the  TEz-case  there  may  be  a  peak  or  a  null 
in  the  broadside  direction  (4>r=0°)  and  the  number  of  lobes  increases 
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with  an  increase  in  the  patch  length.  Smooth  radiation  patterns 
with  single  lobe  can  be  obtained  for  the  TM  -case  with  an  inductive 
surface,  and  the  TE  -case  with  a  capaciti vezsurface.  It  is  also 
interesting  to  point  out  that  there  are  lobes  in  the  deep  shadow 
region  of  the  radiation  pattern.  These  lobes  are  created  by  the 
diffraction  of  the  waves  that  creep  around  the  conducting  portion 
of  the  cylinder,  and  the  number  of  lobes  in  the  deep  shadow  region 
depends  on  the  radius  of  the  cylinder  rather  than  the  impedance 
patch  length.  The  number  of  lobes  in  the  shadow  region  increases 
as  the  radius  of  the  cylinder  becomes  larger;  however,  the  level 
of  the  fields  in  that  shadow  region  decreases  considerably  with 
an  increase  in  the  cylinder  size  to  the  point  where  the  strength 
of  these  lobes  is  negligible  in  comparison  to  the  fields  in  the 
lit  region. 

The  reflection  coefficient,  R  can  be  used  as  a  gauge  for 

the  lobe  sizes,  i.e.,  the  ratio  between  the  lobe  maxima  and  minima. 
For  values  of  f R  |  <  0.9  which  is  the  region  of  interest  in  this 

work,  the  lobe  size  in  the  lit  region  reduces  as  the  magnitude 
of  the  reflection  coefficient  decreases.  In  other  words,  antenna 
structures  with  high  R  leads  to  a  larger  break  up  in  the  radi¬ 
ation  patterns;  this  aspect  will  be  discussed  further  at  the  end 
of  this  chapter. 

Additional  useful  information  about  the  radiation  pattern 
can  be  obtained  if  one  studies  the  behavior  of  the  far  field  pat¬ 
tern  of  the  antenna  at  specific  observation  angles;  namely  4>R=0 
and  <j)n=90  as  a  function  of  the  value  of  the  surface  impedance, 
and  also  separately  as  a  function  of  the  length  of  the  impedance 
patch.  In  Figures  36  through  41  the  magnitude  of  the  field  u(0) 
at  <j)n=0°,  and  likewise  u ( 90 )  at  4>R=90°  are  plotted  in  dB  for  both 
the  TE  and  the  TM  cases.  The  magnitude  of  u(90)  in  dB  for  the 
relevant  unperturbed  geometries  (perfectly  conducting  cylinders 
without  any  impedance  patch)  are  also  included  in  these  plots  for 
comparison.  The  plots  are  normalized  such  that  )u(0)|  for  the 
unperturbed  geometry  is  taken  as  zero  dB. 

A  striking  conclusion  one  draws  from  these  data  is  that, 
for  the  TM  case  with  a  capacitive  impedance  patch,  there  are  very 
sharp  maximums  which  are  spaced  every  half  electric  wavelength, 

X  /2.  In  such  resonances,  one  may  realize  up  to  15  dB  increase 
in  the  magnitude  of  the  field  in  the  end  fire  direction  (<j>„=90°) 
over  that  which  would  exist  for  the  unperturbed  geometry.  *The  TE 
case  with  an  inductive  impedance  patch  also  exhibits  an  increase  z 
in  the  field  strength  in  the  end  fire  direction  compared  to  that 
of  the  unperturbed  geometry.  However,  the  sharp  resonances  of 
the  TM  case  do  not  appear  to  be  present  in  the  TE  case.  Some 
typical  radiation  patterns  for  the  TE  and  TM  cases  are  presented 
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Figure  37--The  magnitude  of  the  radiation  field  of  a  line  source 
(illustrated  in  Figure  10)  as  a  function  of  the  half 
patch  size,  a4>  /X,  . 
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Figure  38--The  magnitude  of  the  radiation  field  of  a  line  dipole 
source  (illustrated  in  Figure  10)  as  a  function  of  the 
half  patch  size,  a4>  /*h> 
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Figure  41— The  magnitude  of  the  radiation  field  of  a  line  dipole 
source  (illustrated  in  Figure  10)  as  a  function  of  the 
normalized  surface  impedance. 


in  Figures  42  through  45  for  both  the  capacitive  and  the  inductive 
type  of  impedances  as  well  as  for  the  case  when  there  is  no  im¬ 
pedance  patch  (unperturbed  geometry).  The  plots  are  normalized 
such  that  their  respective  maximums  are  taken  as  zero  dB. 


Finally,  one  may  point  out  that  the  value  of  the  surface  im¬ 
pedance  may  be  different  or  it  may  be  the  same  for  different  ori¬ 
entations  of  the  sources.  An  example  for  the  first  type  of  surface 
may  be  a  thin  material  of  thickness,  d,  relative  permittivity, 
e  and  relative  permeability,  u  placed  on  a  perfectly  conducting 
surface.  From  the  flat  surfacergeometry,  it  can  be  shown  that [23], 
the  equivalent  surface  impedance  for  the  TE  and  TM  cases  are 
given  as  follows:  z  z 


Zs  =  j  Zo  for  TErcase 


Zs  =  -j  ZQ  for  TMz-case 


(41a) 

(41b) 


where  Vh  is  the  solution  of  the  following  pair  of  equations 


-ud*cotan(ud)  =  urVhd 


(42a) 


(ud)2  +  (Vhd)2  =  (cr-l)(kd)2 


(42b) 


Similarly-  V$  is  obtained  by  solving  the  following  equations  simul¬ 


taneously 

ud*tan(ud) 


e  V  d 
r  s 


(43a) 


(ud)2  +  (Vsd)2  =  (er-l)(kd)2  (43b) 

A  proper  combination  of  d,  e  and  p  may  be  used  to  build  a  surface 
which  can  support  Elliot  type  modes  for  both  the  TE  and  TM  cases 
with  desired  surface  wave  propagation  constants.  Such  a  surface 
may  be  referred  to  as  an  "arbitrary  polarization  surface  wave  struc 
ture"[24].  In  Figure  46,  the  radiation  patterns  of  a  magnetic 
line  source  and  a  magnetic  line  dipole  source  are  presented.  The 
source  is  assumed  to  be  located  on  an  impedance  surface  which  ap¬ 
pears  as  a  capacitive  surface  (Z  =-j710)  for  the  TM  -case  and  an 
inductive  surface  (Z  =j200)  for  the  TM  -case,  so  that  the  surface 
wave  propagation  constants  associated  with  both  cases  are  the  same 
(i.e.,  B=B  =B.=1.13k).  The  radiation  patterns  for  Z  =j710  and 
Z  =-J200  for  the  TM  and  TE  case,  respectively,  aresalso  included 
in  this  figure  for  comparison.  The  plots  are  normalized  to  their 
respective  peak  values. 
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Figure  42--Radiation  pattern  of  a  magnetic  line  source  on  a 
perfectly-conducting  cylinder  which  is  partly  covered 
with  an  impedance  surface  patch. 
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Figure  43--Radiation  pattern  of  a  magnetic  line  source  on  a 
perfectly-conducting  cylinder  which  is  partly  covered 
with  an  impedance  surface  patch. 
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Figure  44--Radiation  pattern  of  a  magnetic  line  dipole  source  on  a 
perfectly-conducting  cylinder  which  is  partly  covered 
with  an  impedance  surface  patch. 


IMPEDANCE 

SURFACE 


SOURCI 


Figure  46--Radiation  pattern  of  a  magnetic  line  source  and 
magnetic  line  dipole  source  on  a  perfectly  conducting 
circular  cylinder  which  is  partly  covered  with  an 
impedance  surface  patch.  The  value  of  the 
impedance  is  assumed  to  depend  on  the 
type  of  the  source. 


In  Figures  47  through  51,  additional  radiation  patterns  of 
a  magnetic  line  dipole  source  (TM  -case)  are  presented  for  values 
of  the  reflection  coefficient,  Rz  moderatly  small.  It  should 
be  noted  that,  these  patterns  aresmuch  smoother  than  the  ones  pre¬ 
sented  in  Chapter  III.  As  indicated  earlier,  the  radiation  pat¬ 
terns  in  Chapter  III  for  the  TM  case  are  included  in  this  work 
only  for  comparison  purposes.  z 

As  indicated  earlier,  the  SCIM  can  be  employed  to  analyze 
the  radiation  pattern  of  a  source  on  a  general  convex  cylinder. 

In  Chapter  II,  the  way  to  extend  the  circular  cylinder  results 
to  an  arbitrary  convex  cylinder  is  presented.  In  Figure  52,  this 
procedure  is  used  in  order  to  obtain  the  radiation  pattern  of  a 
magnetic  line  source  located  on  an  elliptic  cylinder  which  is  covered 
with  an  impedance  surface  patch. 

In  conclusion,  the  present  analysis  of  the  radiation  from 
a  line  source  on  an  impedance  surface  patch  which  partly  covers 
a  perfectly  conducting  convex  cylinder  indicates  that  for  impedance 
surfaces  which  support  an  Elliot  type  mode,  a  moderate  increase 
can  be  achieved  in  the  radiation  intensity  near  the  horizon  (<J>B=90°) 
over  that  which  would  exist  in  the  absence  of  the  impedance  parch. 

In  general,  the  radiation  patterns  associated  with  those  impedance 
surfaces  which  support  the  Elliot  mode  are  found  to  be  rather  sen¬ 
sitive  to  changes  in  the  operating  frequency.  Furthermore,  the 
increase  in  the  radiation  intensity  near  the  horizon  is  accompanied 
by  a  ripple  in  the  radiation  pattern  within  the  lit  region.  The 
addition  of  a  small  loss  and/or  taper  to  the  impedance  patch  would 
decrease  this  ripple  somewhat  without  seriously  degrading  the  radi¬ 
ation  intensity  near  the  horizon.  In  particular,  the  field  strength 
at  the  edges  of  the  impedance  patch  can  be  controlled  by  a  small 
loss,  this  in  turn  allows  one  to  reduce  the  ripple  size.  The  loss 
can  be  introduced  to  the  solution  as  a  perturbation  to  the  loss¬ 
less  case  via  an  exponential  decay  factor  along  the  propagation 
path  on  the  impedance  surface.  It  is  observed  that  the  ripple 
size,  besides  being  dependent  on  other  parameters  (such  as  patch 
length,  radius  of  the  cylinder)  is  also  related  to  the  magnitude 
of  the  reflection  coefficient  [R  |.  As  a  rough  estimate,  an  average 

reduction  of  8-2  dB  in  the  ripple  size  can  be  achieved  upon  re¬ 
ducing  the  magnitude  of  the  reflection  coefficient  by  0.2  in  the 
region  0.5  <  [R  L  <  0.9;  whereas,  it  is  possible  to  smooth  out 
the  ripples  by  4-2  dB  with  a  decrease  in  the  magnitude  of  the  re¬ 
flection  coefficient  by  0.2  in  the  region  ( R,  (  <  0.5.  For  the 
TM  -case,  the  average  ripple  size  at  [R  |  =  0.7  and  |R  |  =  0.9 
isz8  dB  and  18  dB,  respectively.  Similarly,  for  the  TE  -case, 
the  average  ripple  size  at  j RJ  =  o.l  and  |  Rj  =  0.3  is^S  dB  and 
12  dB,  respectively.  It  may  be  possible  to  reduce  the  reflection 
coefficient  by  tapering  the  impedance  patch  without  drastically 
changing  the  actual  antenna  structure,  thereby  reducing  the  ripple 
size  in  the  radiation  pattern  within  the  lit  region. 
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Figure  47--Radiation  pattern  of  a  magnetic  line  dipole  source 
on  a  perfectly-conducting  circular  cylinder  which  is 
partly  covered  with  an  impedance  surface  patch. 
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Figure  49--Radiation  pattern  of  a  magnetic  line  dipole  source 
on  a  perfectly-conducting  circular  cylinder  which  is 
partly  covered  with  an  impedance  surface  patch. 
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Figure  52--Radiation  pattern  of  a  magnetic  line  source  on  a 
perfectly-conducting  elliptic  cylinder  which  is  partly 
covered  with  an  impedance  surface  patch. 


The  present  analysis  based  on  the  SC IM  constitutes  an  exten¬ 
sion  and  improvement  over  the  earlier  GTD  analysis  [2]  which  was 
restricted  to  only  those  impedance  surfaces  which  support  an  Elliot 
type  mode  and  for  which  the  Elliot  mode  propagates  with  negligible 
leakage  on  the  curved  impedance  patch.  These  limitations  of  the 
GTD  solution  are  absent  in  the  SCIM.  While  the  SC IM  in  it's  present 
form  appears  to  be  quite  accurate  in  it's  fairly  broad  regions 
of  validity,  further  work  is  recommended  to  remove  some  of  the 
minor  restrictions  which  currently  exist  (as  outlined  earlier) 
and  also  make  this  solution  applicable  to  a  wider  class  of  problems 
than  the  one  treated  here. 


CHAPTER  V 
CONCLUSIONS 


A  high  frequency  solution  is  developed  for  predicting  the 
radiation  patterns  of  sources  on  a  perfectly  conducting  convex 
cylinder  which  is  partly  covered  by  an  impedance  surface  patch. 

It  is  expected  that  this  study  will  be  useful  to  the  design  of 
on  aircraft  antenna  which  can  be  made  to  radiate  strongly  in  the 
vicinity  of  the  horizon  by  impedance  loading.  Such  an  antenna 
structure  excited  by  a  crossed  slot  on  an  aircraft  fuselage  can 
find  applications  for  satellite  communications.  Thus,  it  is  im¬ 
portant  to  develop  a  theoretical  model  for  this  type  of  an  antenna 
in  order  to  gain  a  better  understanding  of  the  effect  of  the  im¬ 
pedance  loading  on  its  radiation  pattern.  It  is  also  interesting 
to  point  out  that  this  theoretical  model  may  be  useful  in  the  design 
of  antennas  which  will  radiate  strongly  only  in  the  broadside  di¬ 
rection  and  shed  little  energy  in  the  horizon  and  in  the  shadow 
region  of  the  antenna. 

The  surface  impedance  concept  is  an  approximation  for  certain 
types  of  surfaces,  such  as  thin  dielectric  layers,  corrugated  sur¬ 
faces,  slightly  rough  surfaces,  and  imperfectly  conducting  surfaces. 
The  value  of  this  equivalent  surface  impedance  associated  with 
these  surfaces  may  depend  on  the  orientations  of  the  source.  It 
is  pointed  out  earlier  that  the  arbitrary  polarization  surface 
wave  structure  could  be  an  example  for  this  type  of  impedance  sur¬ 
faces. 

The  analysis  of  the  two  dimensional  (2-D)  problem  of  the  radi¬ 
ation  by  a  magnetic  line  source  (TE  case)  and  a  magnetic  line 
dipole  source  (TM  case)  on  a  perfectly  conducting  circular  cylinder 
partly  covered  with  an  impedance  surface  patch  is  employed  in  this 
work  to  simulate  the  3-D  crossed  slot  configuration.  It  is  noted 
that  the  radiation  pattern  of  the  3-D  configuration  would  be  the 
same  as  the  one  for  the  simpler  2-D  case  in  the  roll  plane  con¬ 
taining  the  slots,  provided  the  effect  of  the  wings  are  ignored. 

The  basic  approach  applied  here  for  the  solution  of  the  2-D  prob¬ 
lem  is  called  the  surface  current  integration  method  (SCIM).  The 
SC IM  is  an  important  extension  and  generalization  of  some  of  the 
earlier  solutions  based  on  the  Geometrical  Theory  of  Diffraction 
(GTD)  and  the  hybrid  moment  method  (MM)-GTD  technique  for  analyzing 
this  problem.  The  previous  GTD  solution  was  restricted  to  those 


types  of  impedances  which  support  an  Elliot  type  mode  for  a  given 
source  orientation;  whereas,  the  SCIM  can  be  used  for  an  arbitrary 
uniform  impedance  surface  for  both  the  TE  and  the  TM  cases. 

The  MM-GTD  solution  obtained  previously  fSr  the  TE  case  is  suf¬ 
ficiently  general  (includes  non-uniform  impedance  surfaces)  and 
is  accurate;  however,  the  SCIM  leads  to  a  far  simpler  and  more 
efficient  solution.  On  the  other  hand  the  MM-GTD  solution  for 
the  TM  case  contains  integrals  which  are  very  difficult  to  evaluate 
numerically;  whereas,  in  the  SCIM,  the  TM  case  can  be  handled 
as  easily  as  the  TEz  case.  2 

In  the  SCIM,  one  integrates  the  equivalent  surface  current 
densities  on  the  impedance  patch  in  conjunction  with  an  accurate, 
asymptotic  form  of  the  perfectly  conducting,  convex  cylinder  Green's 
function.  The  relations  between  the  far  field  radiation  pattern 
of  the  sources  and  these  integrals  are  obtained  via  a  generalized 
reaction  theorem  as  explained  in  section  IIA.  The  form  of  the 
equivalent  surface  current  densities  are  found  from  the  canonical 
problems  of  Appendix  B.  In  these  canonical  problems  for  the  TE 
and  the  TM  cases,  it  is  assumed  that  the  cylinder  is  completely 
covered  byzan  impedance  surface.  The  finite  length  of  the  impedance 
patch  is,  then,  introduced  later  into  the  solution  by  a  self-con¬ 
sistent  procedure  for  summing  up  the  fields  reflected  from  the 
edges  of  the  patch  as  presented  in  section  I  IB . 

The  SCIM  can  also  be  used  in  the  solution  of  crossed-slot 
antennas  located  on  an  arbitrary  polarization  surface  wave  struc¬ 
ture  as  long  as  the  coupling  between  the  TE  and  TM  surface  wave 
(or  Elliot)  modes  is  negligible  at  the  edgezof  the  structure  so 
that  the  performance  of  each  of  the  antennas  can  be  analyzed  in¬ 
dependently  as  in  the  ordinary,  single  valued  impedance  surface 
case  and  these  results  can  then  be  superimposed  to  yield  the  total 
radiation  pattern. 

The  effect  of  the  different  parameters  on  the  radiation  pat¬ 
tern  for  both  the  TMz  and  the  TEZ  cases  can  be  summarized  as  follows: 

TMz  Case: 

1)  The  capacitive  impedance  surface  (i.e.,  Z  =- j | x  |  where 
x  is  a  real  number)  supports  an  Elliot  mode  fieldsfor  this  polari¬ 
zation;  hence,  it  is  possible  for  this  case  to  guide  more  energy 
into  the  end  fire  direction  (4>R=90°)  i.e.,  near  the  horizon.  It 
is  observed  for  the  Elliot  mode  case  that  the  antenna  patterns 
become  increasingly  frequency  sensitive  for  those  values  of  Z$ 

=-j|x  |  for  which  |xc|  is  small. 
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2)  The  Elliot  mode  reflection  coefficient,  R  approaches 
unity  as  the  capacitive  impedance  surface  tends  tosthe  perfectly 
conducting  case  (i.e.,  as  |x  |  becomes  small).  Antenna  structures 
with  high  reflection  coefficients  generate  badly  broken  radiation 
patterns  with  deep  ripples. 

3)  An  inductive  impedance  surface  does  not  support  an  Elliot 
mode,  hence  the  field  strength  in  the  vicinity  of  the  horizon  is 
much  lower  than  that  obtained  with  a  capacitive  impedance  surface. 

4)  In  the  case  of  the  capacitive  impedance  patch,  the  number 
of  lobes  within  the  lit  region  increases  as  the  width  of  the  patch 
increases. 

5)  Inductive  impedance  patches  generate  smooth  radiation 
patterns. 

TEz  Case: 

1)  The  inductive  impedance  surface  (Z  =j|x|_|)  supports  an 
Elliot  mode  field  for  this  polarization.  Aft  increase  in  the  field 
strength  near  the  horizon  is  possible  for  this  type  of  impedance. 

2)  The  Elliot  mode  reflection  coefficient,  R.  increases  as 
the  value  of  the  inductive  impedance  increases.  Hence,  inductive 
impedance  surfaces  with  small  |x^|  lead  to  smooth  radiation  patterns. 

3)  The  capacitive  impedance  surface  does  not  support  an  Elliot 
mode,  hence  the  field  strength  in  the  endfire  (horizon)  direction 

is  much  lower  than  that  obtained  with  an  inductive  impedance  surface. 

4)  In  the  case  of  the  inductive  impedance  patch,  the  number 
of  lobes  within  the  lit  region  increases  as  the  width  of  the  patch 
increases. 

5)  Capacitive  impedance  patches  generate  smooth  radiation 
patterns. 
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APPENDIX  A 

AN  ANALYSIS  OF  THE  RADIATION  FROM  LINE  SOURCES  ON  AN 
INFINITE  PLANAR  IMPEDANCE  SURFACE 


In  this  appendix,  the  radiation  from  line  sources  on  an  infi¬ 
nite  planar  impedance  surface  is  analyzed.  The  geometrical  config¬ 
uration  of  the  problem  is  illustrated  in  Figure  53.  Two  different 
source  orientations  are  discussed.  The  first  of  these  is  for  the 
case  of  a  z-directed  magnetic  line  source  excitation,  which  will 
henceforth  be  referred  to  as  the  TE  case;  this  case  is  treated^ 
in  Section  I  below.  Likewise,  the  treatment  of  the  case  of  an  x- 
directed  magnetic  line  dipole  source  type  excitation,  which  is  hence 
forth  designated  as  the  TM  case,  is  dealt  with  subsequently  in 
Section  II. 


Ay  MAGNETIC  LINE  SOURCE 


O' 
(o.o) 
X 


b-r  Ifllllifllitll  !  If  I  I  n  >  rim  rrfn  rrrrrl 


\ 


IMPEDANCE  SURFACE 


(a)  TEz  -  CASE 


(b)  TMz  -  CASE 

Figure  53— Geometry  of  the  canonical  problem  of  a  line  source 
on  an  infinite,  planar  impedance  surface. 
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I.  TEZ  or  Magnetic  Line  Source  Excitation  Case 


The  magnetic  source  density  for  this  polarization  is  denoted 
by  M;  it  is  given  by 

=  2  M  6(x)  5(y)  (A-l) 

where  M  is  the  strength  of  the  source  and  is  assumed  known.  The 
magnetic  field  generated  by  such  a  source  has  only  a  z-component 
and  it  will  be  denoted  by  H  .  H  satisfies  the  reduced,  inhomo¬ 
geneous  wave  equation,  and  the  following  boundary  conditions: 

(V^+k2)Hz  =  jkY0M6(x)6(y)  ;  y>0,  |x|<«  (A-2) 


3Hz 

+  ahHz  =  0  at  y  =  0  (A-3) 

where 

Zs 

ah  =  - jk  f  (A-4) 

o 

For  the  problem  of  interest,  is  chosen  to  be  a  real  number  (i.e., 
Re  Z  =0);  however,  the  method  of  solution  is  applicable  to  complex 
a.  as  well.  The  impedance  boundary  condition  given  in  (A-3)  can 
asymptotically  be  obtained  from  the  more  general  form  of  the  bound¬ 
ary  condition  given  in  (5a)  with  the  aid  of  the  Maxwell's  equations. 
Z  =1/Y  is  the  surface  impedance,  Z  =  Y'1  is  the  free  space  im¬ 
pedance.  is  the  two  dimensional  Laplacjan  operator.  Also,  H 
satisfies  tne  radiation  condition  for  an  eJwt  time  dependence.  1 

In  the  process  of  solving  (A-2)  one  defines  a  two-dimensional 
Green's  function,  G(x,y|x',y')  such  that 

H z(x,y)  =  -jkY0M  Gh(x,y|0,0).  (A-5) 


G  (x,y|x',y')  is  the  solution  of  (A-2)  for  a  -S(x-x' )6(y-y' )  type 
source  term  on  the  right.  It  also  satisfies  (A-3)  and  the  radi¬ 
ation  condition.  Using  the  procedure  for  constructing  higher  dimen¬ 
sional  Green's  function  as  in[25],  one  obtains. 


h  i  °°-Je 

Gh(x,y|0,0)  / 


- jjk2-k2y- jk  I  x  I 


-®-Je 


(ik2-k2  + 


J«h) 


dkx> 


(A-6) 
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where  e  is  a  positive  but  otherwise  arbitrarily  small  number.  The 
proper  kx  sheet  associated  with  the  integrand  of  (A-6)  is  one  for 

which  ImJ<2-kx<0  with  branch  points  at  ±k.  The  integrand  also  has 

poles  located  at  k*  =t/k^+a?  on  the  proper  sheet  for  Im  Zs>0.  These 
properties  of  (A-3r  are  illustrated  in  Figure  54. 


Im  kx 


Figure  54--Contour  of  integration  in  the  complex  kx-plane. 

The  expression  in  (A-6)  can  also  be  represented  as  a  perturbation 
jf  the  fields  on  the  perfectly  conducting  surface  as  follows: 


Gh  =  gJ!  +  gJ 
c  p 


(A-7) 


where  G"  is  the  green's  function  for  the  perfectly  conducting  ground 
plane;  Tt  is  given  by 


,  ,  ». je  -Jkx! xl 

L.<,h—2  dkx  *2j 


’-•te  Jk^i 


(A-8) 


G  is  the  perturbation  Green's  function  and  has  the  following  form 
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! 


G 


h 

P 


( A-9) 


In  both  (A-8)  and  (A-9),  the  observation  point  is  brought  on 
to  the  impedance  surface  (i.e.,  y=0). 


The  two  alternative  expressions  for  the  Green's  function  given 
by  (A-6)  and  (A-7)  may  be  approximated  asymptotically.  The  asymp¬ 
totic  evaluation  of  (A-6)  is  discussed  below  in  part  I-a  and  that 
of  (A-7)  is  discussed  in  part  I-b. 


a.  Saddle  Point  Approximation 

In  the  asymptotic  evaluation  of  (A-6),  it  is  convenient  to 
begin  by  introducing  the  followina  polar  transformations: 


kx  =  k  cosn  ;  dkx  =  -k  sinp  dn  (A-10) 

and 

x  =  pcos<j>  ;  y  =  ps  i  n<t>  (A-ll) 


The  transformation  in  (A-ll)  is  shown  in  Figure  55.  The  Green's 
function  Gn  of  (A-6)  is  now  expressed  in  the  complex  n  plane  as 

h  i  f(n)ft 

G  -  -  - —  /  F.(n)e  dn  (A-12) 

2itj  c 

where 

Fh(n)  =  - rr  (a-13) 

sinp+ja^k 

and 

f(n)  =  -j  cos(n-4>) ;  0=  kp.  (A-14a); 

(A-14b) 

fi  is  the  large  parameter  in  this  asymptotic  development.  The  contour 
of  integration  c  is  indicated  in  Figure  56. 

nD  and  nD  are  the  poles  of  (A-6)  in  the  complex  n  plane. 

*1  k2 

By  setting  the  denominator  of  (A-13)  to  zero,  one  can  easily  show 


86 


Figure  55--Polar  and  rectangular  coordinate  systems 


Im  r) 


Figure  56--Contour  of  integration  in  the  complex  n -plane 


that  r^,  =-j  sinh’^(ah/k)  and  np  =7T+j  sinh"^  (ci^/k ) .  The  saddle 

point  is  denoted  by  n,  and  it  is  the  solution  of  f'(n)=0.  It 
can  be  shown  that  n  =  4>  for  (0<4><tt).  For  the  special  case  when 
the  observation  point  is  on  the  impedance  surface  (i.e.,  4>=0), 
ns=0. 

Using  the  results  of  Felsen  and  Marcuvitz[26]  one  can  obtain 
an  asymptotic  solution  for  (A-6)  as  follows 


Gh  ~  Re  e' 


pf(n) 


2tt  j 


±j2aJ7e‘Qb2  Q(+jbJfi)  +  J|  | 


( A- 1 5 ) 


Im  b  <  0, 


where 


R  =  1  im 

rp-nr 


(n-np^)Fh(n)e 


-f(n)n| 


,  a  =  1 im  [(n-np  )Fh(n)]  (A-15a) 
rr*nPi  1 


T(n! 


f(n  )  ,  e 
P1 


fll  ,  Jim  b> 0^* 
Of  for  I  Im  b<0 


(A-15b) 


Q(y)  =  /  e"x  dx 
y 


Y~  erfc(y)  (A-15c) 

where  erfc  is  the  complementary  error  function.  Alternatively,  sub¬ 
stituting  (A-13)  and  (A-14)  into  (A-15),  with  ns=0,  yields 

2 

lh  _Jphx  ikx  I  +  ^h  -  i—  _Kxr 

_ _!I  a  n  4.  o  I  I  _ [L  0^4-tK 

6, 


+  -  i—  _kxbh 

"  Q(+Jbh^)e 


1 


TT  ^h 


2rr  j  J  kx  &  ,b 


h  hJ 


Im  b^<  0 


(A-16) 


where 


h  .  /7  ah  /  Bh  ,  n-*« 

bh  -  r  (  r  +  1} 


h2  * 


(A-16a) 
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for  (a^<0)  or  ( Im  b^<0). 


( A-  16b ) 


b.  Asymptotic  Evaluation  of  the  Perturbation  Green's 
Function  Representation 


The  first  term  in  (A-7)  was  readily  evaluated  in  (A-8)  as  the 
Green's  function  for  the  half-space  with  a  magnetic  line  source 
located  at  the  coordinate  origin.  The  effect  of  the  impedance  is 
represented  by  the  perturbation  Green's  function,  G"  which  is  given 
in  (A-9).  As  indicated  earlier,  the  integrand  in  (fl-9)  has  poles 
for  positive  values  of  ot,.  In  the  process  of  evaluating  (A-9), 
one  defines  a  new  quantity  as  follows 


-jo 


-  1 
h  2ttj 


where  ah  =  "  lahl  • 

Through  a  simple  manipulation,  (A-9)  becomes 


-2jo 


h  Wj 


! 


-je 


-j*xM 


'-“-je  (k2+a^)-kx 


2  dkx  +  Gp 


( A- 1 7 ) 


( A- 18 ) 


The  first  term  in  RHS  of  (A-18)  can  be  calculated  by  enclosing 
the  contour  in  the  lower  half  complex  k  -plane  and  evaluating  the 
residues  enclosed.  Hence,  (A-18)  becomes, 


G 


h 

P 


(A-19) 


The  integrand  in  G 
it  can  be  rewritten  as 


!?  has  no  poles  on  the  proper  sheet,  therefore 
follows: 


J  H^2)(kjx2+x2)  e“hT  dr  (A-20) 

where  the  relations 

1  *  -jdk2-k2  +  ja"3r 

-= -  =  j  J  e  dx 

^h 
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and 


-jkx  |x | - j T 


k2-k2 

X 


are  employed  [26]  in  the  integrand  of  (A-17)  to  arrive  at  (A-20). 

For  large  values  of  aT  and  x,  (A-20)  can  be  evaluated  by  the 
large  argument  approximation  for  the  Hankel  function  as 


,h\  ah 


2j 


f  0 


-j(kx  [l+  j(  1  )2]  -  +  a^x 


uk  jx  ] 


V 


dT.  (A-21 ) 


Next,  introducing  a  new  variable  z  =  Tjk" (  ~  +  j rr^  U  \  one  can 
chnw  that  v  PT*  k  -J  2  J 


show  that 


JYt 


oj  t  -  n  c 

P  •  —3.  e 

kJiT 


2 

fa  -Jk|x|(i  -  ~2 


|x|-.  (A-22) 


For  small  values  of  a^,  the  integral  in  (A-20)  is  slowly  conver¬ 
gent;  however,  one  could  overcome  this  difficulty  by  expanding  Gh  in 
a  power  series  of  a”  as  follows.  p 


1 

2tt  j 


j-je 

-®-je 


dkx 


x| 


__1 _ 

(k2+a2)-k2 
rr  x 


H _ 


(A-23) 


Evaluating  the  first  term  in  (A-23)  by  the  calculus  of  residues,  and 
expanding  the  second  term  in  powers  of  ct^  yields 
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-JBjxl 


.  ah  "Jph 

•Ja:e 


■je  -Jky|x| 


h  uh  f  dk  e  M"'x 1 

-  ?ij  *co.je  dkxe 


1  ah  ,  ah 

raw  (kz.k2)V2 


+  . . . 


(*C-WC  (k‘-k‘) 


(A-24) 


or 


,h- 


‘J*h 


rmU^  H> (k|x|) 

( A-25 ) 


In  obtaining  (A-25)  from  (A-24),  the  relation 

\P 


st?  j e  2  2  TFr^ dk,t  =  ^  T^^rferT  Hp2)<kW:  p=1,2,3'" 

--*(k2-k2)P.P^  (A-26) 


is  used. 


A  closed  form  expression  for  (A-25)  may  be  obtained  for  large 
values  of  x.  To  this  end,  one  uses  the  asymptotic  approximation 
for  the  Hankel  function  and  rearranges  the  resulting  series  in  (A- 
25)  which  yields 


e 


-jBhlxl  j<*h 
+  W~ 


-  jk  |  x 


e 


J  a\  k|*| 
.  2k‘ 


or 


2n 

1  3. . . (2n+l) 


(A-27a) 


(0H 


(A-27b) 
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Since  (A-27b)  is  derived  for  small  values  of  a. ,  one  can  approxi¬ 
mate  Sh  -  k  for  the  ampl ilude-terms  (i.e.,  when  0.  ts  in  the  denomi¬ 
nator),  and  by  3h  ~  k(l+ar/2k^)  for  the  phase  terms  (i.e.,  when 
it  is  in  the  exponent).  This  approximation  yields 


(A-28) 


.-  The  above  analysis  shows  that  regardless  of  the  value  of  ah, 
G"  has  the  same  form  as  given  by  (A-22)  or  (A-28),  for  large  val 
of  | x  |  • 

Substituting  (A-27)  or  (A-28)  into  (A-19)  and  substituting 
this  result  and  (A-8)  into  (A-7)  yields 


ues 


G  ~  -  a —  e 
Ph 


£h  + 


ja. 


k  JiT 


,-jk|x| 


frrrrV-klxlBh 


|(j  Bjkfxj^e 


+  h 


M) 


(k|x|) 


(A-29) 


where 


Bh 


(A-29a) 


As  indicated  earlier,  (A-29)  is  valid  for  both  small  and  large 
values  of  a.,  when  |x|  is  large.  It  can  easily  be  shown  that  (A- 
29)  goes  tonthe  saddle  point  approximation  of  (A-16)  for  small  values 
of  a.  and  large  values  of  |x|.  On  the  other  hand,  for  small  values 
of  | x | .  the  Hankel  function  in  (A-29)  will  dominate  and  provide 
an  adequate  description  of  the  current  near  the  source.  Thus  Eq¬ 
uation  (A-29)  is  an  asymptotic  approximation  of  the  Green's  function 
which  also  remains  valid  in  the  vicinity  of  the  source  for  all  values 
of  a. ;  whereas,  (A-16)  is  an  asymptotic  approximation  of  the  Green's 
function  for  large  values  of  x. 

II.  TMz  or  Magnetic  Line  Dipole  Source  Excitation  Case 


The  magnetic  source  density  for  this  polarization  is  denoted 
by  ffld,  and  it  is  given  by 


where  Md  Is  the  strength  of  the  source  which  is  assumed  known. 

The  tangential  magnetic  field  generated  by  such  a  source  has  only 
x-component  and  it  will  be  denoted  by  H  .  H  can  be  obtained  from 
the  only  electric  field  E*  generated  byxan  electric  line  source 
I  =  z  IS(x-x')6ly-y').  Here,  I  is  the  known  strength  of  the  electric 
line  source.  E*  satisfies  (A-2)  and  (A-3)  with  the  exception  that 
YqM  Is  replacedzby  21  and  is  replaced  by  a$  respectively,  where 


E*  also  satisfies  the  radiation  condition  for  an  e^a,t  time  depend¬ 
ence.  With  above  changes,  (A-5)  becomes 


E*(x',y')  *  -jkZ0IGs(x',y'|x,y) 


(A-32) 


Gs  corresponds  to  Gh  of  the  TE  -case,  and  using  the  procedure  of  [25], 
one  can  show  that  z 


Gs(x',y'|x,y)  = 


-jk  |x'-x[ 

e  x  dkx;  y'>y  (A-33) 


where  R$  »  (j k2-k*  -  jas)(>/k2-k^  +  ja$)_1. 


The  case  for  y>y'  can  be  obtained  from  (A-33)  by  Interchanging  y 
and  y'.  It  is  no  surprise  that  (A-6)  is  the  special  case  of  lA- 
33)  when  y-0  and  a$  is  replaced  by  c^.  The  magnetic  field,  M^x 

generated  by  T  can  be  found  from  VxE*z*-ju>y  TF1.  The  electric  field, 
E  z  generated  by  Wd  ofjA-30)  can  now  be  found  from  the  reciprocity 
tneorem  which  yields  I  z»E  *  -Mdx*R*.  Making  use  of  the  other 

Maxwell's  curl  equation,  H  *  -( jwygr^VxE+RJ  it  can  be  shown  that 
Hx(x,y)  *  -jkY0Md  Gs(x,y|0,0)  (A-34) 


Gs(x,y|0,0) 


Jas  I 

/ 


Je  ^4 

-jE  JkV  +Ja, 


where 


The  contour  of  integration  in  (A-35)  and  the  branch  cuts  in 
the  intearand  (A-35)  are  the  samj  as  in  Eigure  54.  The  Integrand 
in  (A-35)  has  poles  located  at  k~  =  Jlr-aF  on  the  proper  sheet 
for  the  positive  values  of  a$.  p  s 

The  expression  in  (A-35)  can  also  be  represented  as  a  pertur¬ 
bation  of  the  fields  on  the  perfectly  conducting  plane  surface  as 
follows: 

Gs  =  G*  +  Gp  (A-36) 


where  G*  is  the  TM  type  Green's  function  for  the  perfectly  con¬ 
duct  ingcground  plane;  it  is  given  by 


1 

?ir  j 


i  °°-Je  f 

K  !  J 


-oo- je 


(kj x| ) 


(A— 37) 


Gs  is  the  perturbation  Green's  function,  and  it  has  the  following 
fBrm: 


G 


s 

P 


1 

"5rrj 


'"Je  k2(Jk2-k^+ja$) 


(A-38) 


In  both  (A-37)  and  (A-38),  the  observation  point  is  brought  on  to 
the  impedance  surface  (i.e.,  y=*0).  The  integral  in  (A-35)  will 
next  be  asymptotically  approximated  via  the  saddle  point  method 
in  part  II-(a);  likewise,  the  perturbation  integral  corresponding 
to  Gp  in  (A-38)  will  be  asymptotically  approximated  in  part  II-(b). 


a.  Saddle  Point  Approximation 


Equation  (A-35)  will  be  asymptotically  evaluated  in  a  manner 
similar  to  that  in  part  I-a  of  this  appendix. 


Substituting  the  transformations  (A-10)  and  (A-ll)  into  (A- 
35)  yields 

s  1  ,  f(n)n 

Gs  -  -  ;  Fs(n)  e  dn  (A-39) 


where  .  ? 

Fste)  - 

s  k  slnn+j  / 

and  f(n)  is  the  same  as  in  (A-14). 


(A-40) 
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Figures  55  and  56  are  also  the  same  for  TM  case  with  the_ excep¬ 
tion  that  the  poles  of  (A-40)  are  now  located  at  nD  =  -j*s1nh-1(a./k) 

-1  V1  5 
and  nD  s  *+J  slnh  (ac/k).  Changing  the  subscript  or  superscript 

v2  s 

"h"  to  "s“  in(A-15)  through  (A-15c),  one  arrives  at  an  asymptotic 
representation  of  G.  The  final  result  has  the  following  form 


where 


jotg  -Iklxl  ^°S  -  r . 

*  sj  :  Im  bs?  0 


>r. 


•{!) 


for  (as<0)  or  (Im  bs<0) 


-jkxb. 


(A— 41) 


(A-42a) 


(A-42b) 


b.  Asymptotic  Evaluation  of  the  Perturbation 
6reen‘s  Function  Representation 


The  alternative  representation  in  (A-36)  will  now  be  evaluated. 
It  is  already  shown  In  (A-37)  that  G"!  is  the  Green's  function  for 
the  perfectly  conducting  plane  surface  with  an  x-directed  magnetic 
line  dipole  source  located  at  the  coordinate  region.  The  Impedance 
effect  Is  then  represented  by  G"*  which  Is  given  in  (A-38).  In  the 
process  of  evaluating  (A-38),  oRe  first  separates  the  pole  contri¬ 
bution  as  follows. 


^  L 


k2-k2 
*  Kx 


-jk  J x| 


k2(Jk2-k2+joT) 


»-je  j  I"  2jaJ_ 
-~-je  ?  k2ta?-k^ 


1  -jk  |x| 

-  2 jet  I  e  x  dk 


(A— 43) 


where  a"  =  -|a  | .  The  last  Integral  In  (A-43)  can  be  Integrated 
via  the  calculas  of  residues. 


1 


Next,  Incorporating  the  following  relationships 

j2 


*  Vtmt) 


-jkx|x|  k2-k2  -jkx|x| 
e  =  — s —  e 


and 


k‘ 


Jk2-k^+jax  Jk2-k2  kjk2-k2  [Jk2-k2+jas] 


(A-44a) 


(A-44b) 


into  the  first  integral  of  (A-43)  yields 
jZ 


V  dfklxllv 


-jkx|x| 


h 


-je 


dk. 


*  T  «<klxP=s 


ja|  -J0S  |x  | 
- 2  e  H 

B$kZ 


(A-45) 


Comparison  of  (A-17)  with  the  term  in  the  curly  bracket  of 
(A-45)  shows  that  they  are  the  same  with  the  exception  that  in  the 
latter  cl  is  replaced  by  cl.  Hence  the  evaluation  of  (A-17)  leading 
to  the  result  in  (A-22)  iPequally  valid  for  the  integral  in  (A- 
45)  as  long  as  all  of  the  superscripts  and  subscripts  in  part  I- 
b  are  replaced  by  "s"  in  the  present  development.  With  the  above 
comments  in  mind,  one  can  now  readily  write  an  asymptotic  expression 
for  Gp  by  substituting  (A-22)  or  (A-28)  into  (A-45). 


~(i  +  — — n£2)(k|x|) 

V  d[k|x|]27[2J  0 


-  £^-jk|x|  Q(jb  Jk]7]) 
k  Jtt  5 

2cl 

+  5(k  |x|)e$ 


~k  |x  |b2 


jal  “j0s|x| 

— |  e  «. 

J  y2 


(A-46) 


where 


(A-46a) 
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The  term  Involving  the  Q  function  In  (A-46)  can  be  written  as  follows 


£^s  -jk|x| 

k  yr 


Q(jbjk]7|)e 


Q(z) 


(A-47) 


where  z  *  jb  Jk|x|.  The  term  Inside  the  brackets  Is  recognized 
as  the  large5argument  approximation  of  the  Hankel  function.  Hence 
one  can  conjecture  a  different  asymptotic  form  for  (A-46)  as  follows 


G 


s 

P 


Ho2)<kM) 


1 


6  ( k  |  x  | )  es . 


z2 

(k |x| )2z  ez  Q(z) 


(A-48) 


This  conjecture  can  be  justified  if  one  notes  that  as  the  admit¬ 
tance  tends  to  Infinity  corresponding  to  the  case  of  a  perfectly 

conducting^ground  plane  then  2z  ez  Q(z),  tends  to  unity  (i.e., 

llg^  2z  ez  Q(z)  *  1)  for  xj*0.  Hence,  the  first  two  terms  In  (A- 
s 

48)  cancel  each  other  as  <*•**.  The  addition  of  a  small  loss  in 
the  medium  for  (y  >  0,  |x|<»7  nullifies  the  surface  wave  as  «■*»>, 
thereby  resulting  in  a  zero  contribution  from  G*  term  to  G.  s 
Hence,  In  the  limit  as  the  Impedance  surface  goes  to  the  perfectly 
conducting  ground  plane,  Gs  becomes: 


11m 

0t 

|x|« 


Gs 


,s  _  1  «S2)0<|x| 

c  ^  k  J  x 


(A-49) 


which  Is  the  desired  result  for  TMZ  case. 

Finally  substituting  (A-37)  and  (A-48)  Into  (A-36)  yields: 
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Gs  - 


j°C  -jbsM 

-r-  e  5 


k  p, 


2a  i  Hj^(k|x| ) 

es  +  nr  6(k*x*)cs +  75  TTqiq - 

j2 


H»2,(k|x„  •  f  ■ 32 

1  '  ,J  (A-50) 


where  as  before  z  =  jb/k | x|  and  b  is  defined  in  (A-46a).  The  first 
term  in  (A-50)  corresponds  to  the  Surface  wave  field. 


APPENDIX  B 

AN  ANALYSIS  OF  THE  RADIATION  FROM  LINE  SOURCES  ON  A  CIRCULAR 
CYLINDER  COMPLETELY  COVERED  WITH  AN  IMPEDANCE  SURFACE 


In  this  appendix,  the  radiation  from  line  sources  on  a  circular 
cylinder  covered  with  an  impedance  surface  is  analyzed.  The  geo¬ 
metrical  configuration  of  the  problem  is  illustrated  in  Figure  4. 

Two  different  source  orientations  are  discussed.  The  first 
of  these  is  for  the  case  of  a  z-directed  magnetic  line  source  ex¬ 
citation  which  is  illustrated  in  Figure  4;  this  case  will  be  re¬ 
ferred  to  as  the  TE  case  and  it  is  treated  first  in  Section  I  below. 
The  second  source  orientation  refers  to  a  $' -directed  magnetic  line 
dipole  source  excitation  case  as  shown  in  Figure  4.  Henceforth, 
the  second  case  will  be  designated  as  the  TM  case  and  it  will  be 
dealt  with  subsequently  in  Section  II.  z 

I.  TEz  or  Magnetic  Line  Source  Excitation  Case 


The  magnetic  source  density  for  this  polarization  is  denoted 
by  TT  and  is  given  by 


W  =  z  M  6lP--.a)M<t>l .  (B-l) 

where  M  is  the  known  strength  of  the  source  and  a  is  the  radius 
of  the  cylinder.  The  magnetic  field  generated  by  such  a  source 
has  only  a  z-component  and  it  will  be  denoted  by  H  .  H  satisfies 
the  reduced,  inhomogeneous  wave  equation,  and  the  following  boundary 
conditions: 


(v|  +  k2)Hz  =  jkY0  M  ; 


p>a 

0<4><2tt 


(B-2) 


and 


9H 

sT  +  »  0  »t  p  ■  a  (B-3) 

where  CL^is  given  in  (A-4).  Also  H  satisfies  the  radiation  condition 
for  an  iJU)t  time  dependence.  For  solving  (B-2)  one  defines  a  two- 
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h 

dimensional  Green's  function,  5  ( p,<j> | p '  <j)')  such  that 
Hz(p,4>)  =  -jkYQM  Sh ( p ,4, |  a,0 ) 


(B-4) 


Gh(p Ip  ‘  »4> 1 )  is  the  solution  of  (B-2)  with  -6 (p-p  ’ ) 6 ( -<j>  ’  )/p 
source  term  on  the  RHS.  It  also  satisfies  (B-3)  and  the  ra¬ 


the  source  term  on  the  RHS.  It  also  satisfies  (B-3)  and  the  ra¬ 
diation  condition.  It  can  be  shown  that  the  Green's  function  has 
the  following  form  [26]: 


Gh(p  »4> | P  ‘  *4> ' )  =  -  4  /  ^  du  H(1)(kp  ' )  - 

a  v 


*  Hj2)(kp) 

where  Qh  =  TpTpJ  +  r  • 
One  may  write 


cos  y(  7T-  M 
sin  vit 


(kp'j 


(B-5) 


sinvir 


I  [e‘ 
1=0 


1-jv(2Tr£) 


(B-6) 


where  <|>“=2tt- [(J)-4> '  | .  Physically,  the  above  series  corresponds  to 
multiple  encirclements  of  the  field  around  the  cylinder  in  the  azi¬ 
muthal  propagation  representation  for  Green's  function  of  (B-5). 
Since  one  is  interested  in  applying  the  results  of  the  present  anal¬ 
ysis  to  perfectly  conducting  cylinders  with  an  impedance  surface 
patch,  the  effect  of  these  multiple  encirclements  will  be  neglected 
with  the  exception  of  the  ray  which  will  travel  the  shortest  dis¬ 
tance  between  the  source  and  the  observation  points,  i.e., 

(cosv(it-|<J)-<j)  '|  )/sir\m  will  be  replaced  by  je~Jvl^“^  I  in  (B-5), 
yielding 


G  (p  ,4>  |  p '  ,<|> ' ) 


Je 

;  . 

-oo-Je 


H^2)(kp') 


fl)  (ka) 

dv  Hy>(kP')-  — 

Lv  QhHi  (ka) 


(B-7) 


From  here  on,  Gn  in  (B-7)  will  be  referred  to  as  the  Green's  func¬ 
tion  and  it  will  be  used  in  place  of  (for  example  Equation  (B- 
4)  now  will  read  H  (p,*)  =  - jkY  m  Gh(P  ,<H  a,0)).  One  notes  that 
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the  contributions  from  higher  order  encirclements  to  G  in  (B-5) 
can  be  found  simply  by  replacing  '|  in  G  of  (B-7)  with  1 4> -4> ' | 

+2^  or  with  (2jr- 1 <j>-d>*  | )  +2 ttH  and  summing  them  up  over  l  as  indicated 
by  (8-6). 


The  expression  in  (B-7)  can  be  represented  as  a  perturbation 
of  the  fields  on  the  perfectly  conducting  cylindrical  surface  as 
follows: 


Gh  =  gJ  +  Gj  (B-8) 
c  p 

where  G^1  is  the  Green's  function  for  the  magnetic  line  source  in 
(B-l)  when  it  is  placed  on  a  perfectly  conducting  cylinder;  it  is 
obtained  from  (B-7)  upon  setting  =  0.  Thus, 


°°-je 

Lu <■ 


(d) 


H^2)(ka) 


a-jv|4>| 


(ka) 


(B-9) 


Following  an  analysis  which  is  based  on  some  earlier  work  by 
Hasserjian  and  Ishimaru [19] ,  one  obtains  the  following  approximation 
for  (B-9): 


G|!-^H(2,(ka|*|) 


r3/2 


7 

TO 


C3  + 


where 


(B-10) 


fol lowi: 


is  the  perturbation  Green's  function  and  it  is  given  as 


rh  joth 
Gp  ="TT 


Pv(ka»h)  Pv(ka;0)  e^ltl 


(B-ll) 


where 


P 

v 


(ka;ah) 


H(,2)(ka) 


,(2)' 


(ka)  + 


ah 

r 


H^(ka) 

v 


( B-12 ) 
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The  prime  over  the  Hankel  function  in  (B-12)  indicates  differen¬ 
tiation  with  respect  to  the  argument.  In  both  (B-9)  and  (B-ll), 
the  observation  point  is  brought  on  to  the  surface  of  the  cylinder 
(i.e.,  p=a) . 


In  the  process  of  evaluating  G"  in  (B-ll),  one  encounters  the 
problem  of  expanding  the  Hankel  function  in  an  asymptotic  series 
for  the  different  ranges  of  the  parameters  involved.  n  the  im¬ 
mediate  neighborhood  of  the  source;  that  is,  when  ka | <f>  is  small, 
the  main  contribution  to  the  integral  in  (B-ll)  comes  from  the  - 
values  of  v  which  satisifes  the  inequality  |v-ka|>0(m)  where  2rrr(  = 
ka)  is  taken  as  the  large  parameter.  The  Debye  asymptotic  expansion 
is  the  most  suitable  representation  for  the  Hankel  function  for 
small  values  of  ka|<t>|.  If  k a | <t> |  is  moderately  small  then  one  em¬ 
ploys  the  Watson  approximation  for  the  Hankel  function.  The  Watson 
approximation  corresponds  to  the  values  of  v  such  that  |v-ka|<0(m). 


A  third  representation  referred  to  as  the  creeping  wave  rep¬ 
resentation  for  the  fields  can  be  found  by  evaluating  (B-8)  via 
the  calculus  of  residues.  This  creeping  wave  solution  is  rapidly 
convergent  for  large  values  of  ka 1 | .  One  points  out  that  the  Debye 
approximation  presented  below  goes  to  the  Watson  approximation  for 
v-O(ka),  and  that  the  results  obtained  from  the  Watson  approximation 
for  the  Green's  function  blend  with  the  creeping  wave  representation 
for  moderately  large  values  of  ka | <J> | . 

Equation  (8-11)  will  be  evaluated  firstly  by  the  method  of 
Debye  approximation  in  part  la  of  this  appendix.  The  Watson  approxi¬ 
mation  and  the  creeping  wave  representations  of  (B-ll)  will  be 
carried  out  subsequently  in  parts  (lb)  and  (Ic)  respecti vely. 


a.  Debye  Approximation 


This  asymptotic  expansion  is  valid  for  large  values  of  the 
argument,  ka,  and  the  order,  v,  of  the  Hankel  function  such  that 
|v-ka|>0(m).  Then,  it  can  be  shown  that  [27,28] 


(ka)  - 


-jv(tann-n) 


vtan  n 

L 


(B-13a) 


and 


(ka) 


i  IJKTn  +J  T 

J  irn 


1 


V 


+  0 


( B-l 3b) 
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where  n  is  a  fixed  positive  number  and  cos  n  a  v/ka.  Also,  u. 
and  Vj  are  given  as  1 

ul  =  “  (3  cot  ’l  +  5  cot3n)  ;  vi  =  (9  cot  n  +  7  cot3n) 

( B— 14) 


After  substituting  (B-13)  into  (B-12)  and  rearranging  the  terms,  one 
arrives  as 


pv(ka;ah)~ 


_J _ 

Otu 

sinn+j  y~ 


v^sinri+jct^k'^u^ 

sinn+jahk"^ 


(B-15a) 


Carrying  out  the  division  in  (B-15)  and  retaining  only  the  terms 
up  to  and  including  terms  of  order  1/v  with  respect  to  unity  and 
using  (B-14)  yields. 


Pv(ka:“h>- 


sinn+ja^k 


1  + 


2sin2n 


1 

sinn+j^k"1 


cosn  + 
v 


(B-15b) 


Substituting  (B-15b)  into  (B-ll)  and  performing  the  indicated  multi 
plication  between  P  (ka;ah)  and  Pv(ka;0)  leads  to: 


where 


-Je 
/ 


7k 2 - C2  [i k2 - C2+ J' oth] 


(B-16) 


(B-17) 


and 
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rh  '“h  1 
Gp  j'TSj 


oo-JC 

I  . 


df. 


2\v!--C  +  jah 

{k2-^)2[ik2-^2+jah] 


(B-18) 


A  change  of  variable,  (i.e.,  £  =  v/a)  has  been  introduced  in 
the  above  equations,  and  it  is  recongized  that  cos  n  =  v/ka  and 

sing  =  17:  f}/k.  Comparison  of  (B-17)  and  (A-9)  reveals  that  these 
two  equations  are  indeed  identical  with  the  exception  that  the 
distance,  )xj  between  the  source  and  the  observation  point  on  the 
flat  surface  in  (A-9)  is  now  being  replaced  by  the  arc  length  a|4>| 
on  the  surface  of  the  cylinder.  Hence,  the  same  argument  presented 
in  Appendix  1  can  also  be  used  here  for  (B-17);  this  results  in 


j«h  -jVI*!  jah 


k/n 


i'jka M  Q(jbhJkaj 


)e 


-ka|<|,|b2 


(B-19) 


where,  as  before 
6 


■  Jk2+ah  » 


-  OIL 


(B-19a; 


and 


S  'io  for  v° 


(B-19b ; 


In  order  to  evaluate  G"  ,  one  rewrites  (B-18)  as  follows: 

Pj 


o-JC 


Pi 


-  *  fe  /  ■  *  e 

-o°-  je 


(kV) 


i 


(> 


1  _ 

ik2-^2[Jk2-C2+jah]  - 


(B-20) 


The  first  two  terms  in  (B-20)  can  be  evaluated  using  calculus  of 
residues.  The  last  term  in  (B-20)  can  be  written  in  terms  of  G" 
of  (B-17),  leading  to  po 


(B-21) 


The  explicit  form  of  the  derivative  term  in  (B-2)  can  be  found  by 
simply  differentiating  (B-19)  with  respect  to  a^,  resulting  in 


In  the 

fkalTj) 

i\  ’  u  j  n 

definition  of  the  Q  function  (see  (A-15c)),  and  ^  erfc(z)  *  -  e 
It  may  also  be  pointed  out  that  a^/ah  =  +1;  ah<0. 

b.  Watson  Approximation 

The  Watson  approximation  is  valid  for  the  large  values  of  the 
argument,  ka  and  the  order,  v  of  the  Hankel  functions  such  that  |v- 
ka|  <0(m).  For  such  values  of  v,  it  can  be  shown  that  [28] 


process  of  deriving  (B-22)  one  uses  the  fact  that  Q(jbh 
=  -  4  e  h  which  is  the  consequence  of  the 


Hl2^(ka)~J—  w,(o) 


(B-23c 


_  v-ka 
m 


caV 


■( 


,1/3 


rj 


(B-23c) 


w2(o)  =  /i[Bi ( a)  -j  Ai(a)]  ,  w£(o)  =/n[Bi,(a)  -  j  Ai'(o)] 

(B-23d) 

Ai(o)  and  Bi(o)  are  the  Miller  type  Airy  functions.  Substituting 
(B-23a,b)  into  (B-12)  yields 


%<“*•<*>-  ’T - ~ 

i  Vo)-  r 


(B-24) 


where  <Mo)  the  ratio  of  wo(0)  to  w?(a)  and  has  the  following 
asymptotic  form 

Wo(o) 


1  /i  \ 

Vo)  -*o  +0(^J 


(B-25) 


Substituting  (B-25)  into  (B-24),  expanding  (B-24),  and  retaining 
only  the  terms  of  order  1/ka  leads  to 


P  (kajaJ-  - 


1 


m_1Va-  L 


1  +  7 


(nf1,/a)2[nf1,/o-  f1 


(B-26) 


Substituting  (B-26)  into  (B-ll)  and  keeping  only  the  terms  up  to 
the  order  of  ( — — sA  leads  to 


Gh  ~  Gh  +  Gh  fi~)  +  of— — 

P  P0  PiVkV  V[ka]V 


(B-27) 
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where 


and 


H  e-jka!^l 

-W  ~Td - 


jco 

* 

-Job 


?  a  |4>!S 

e _ 


\  jjj  o^fk 


e-jka  |4»| 

TiT 


ds 


2/s  +  jyjahk_5s 
s2  (/s+j  /jothk_ls)2 


(B-28) 


7al<t»ls 

e 

(B-29) 


In  the  process  of 
change  of  variable  s  = 

a  a  s*  and  natura11y 

Equations  (B-28)  and  (B-29)  are  in  the  form  of  an  inverse 
Laplace  transformation.  Equation  (B-28)  can  be  readily  evaluated 
using  standard  Laplace  transformation  tables  [28j: 


obtaining  (B-28)  and  (B-29)  one  employs  a 
-j2(^  -  k).  In  the  s  domain,  (B-23c)  becomes 

nf  Vo  =  Jj  s/k . 


erfc  (j  e-^HICViV^2 

(3-30) 

r  faz  l  . — 

Since  erfc(-z)  =  2-erfc(z)  and  z  *  j  - 1  L/ka|<|>|,  one  can  re¬ 
write  (B-30)  as  follows  L  J2  k  J 


h  n 
!P„  '  ‘  T 


1a7 


Again,  In  (B-30b)  b 


and 


Comparing  (B-19)  with  (B-30b)  reveals  that  they  both  are  the 
same  with  the  exception  that  the  surface  wave  term  In  (B-19)  Is 
more  accurate  than  the  one  In  (B-30b).  It  can  be  shown  that  (B- 
30b)  Is  the  limiting  case  of  (B-19)  for  small  impedances  (l.e. , 
is  small,  8^  -  k).  The  second  term,  G|J  In  (B-27)  can  easily 

be  calculated  from  (B-29),  If  one  recognizes  that  (B-29)  can  be 
rewritten  as  follows: 
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7Z3 


J® 

J  ds  e' 

-j® 


1  al4>ls| 


^s(/s+j7Jahk-’) 


d 


(B-31) 


Equation  (B-31)  can  now  be  evaluated  with  the  aid  of  a  Laplace 
transformation  table,  yielding 


W;" 

po 


(B-32) 


In  obtaining  (B-32)  from  (B-31),  the  definition  of  G|J  in  (B-30) 

Is  also  used.  The  explicit  form  of  the  derivative  1n°(B-21)  can 
be  found  by  simply  differentiating  (B-30b)  with  respect  to  a. : 

r.  _h  *1  " 


3ou! 


.  in  f  Jfc«UI  ,-Jk*l*l  ;  H*1*1  -jk.|*| 

, _  -JkaUlB? 

<HJSjkjj*|)e  n  ;  ah  S  0  (B-33) 


Comparison  of  (B-21)  and  (B-22)  with  (B-32)  and  (B-33)  reveals 
that  the  latter  pair  Is  the  special  case  of  the  former  pair  for 
small  Impedances,  and  larger  arc  lengths,  a|$|.  It  can  be  seen 
that  In  the  limit  as  a|$|  becomes  large  and  a.  becomes  small,  which 
Is  the  region  where  the  Watson  approximation  fs  valid,  the  Debye 
approximation  goes  to  the  Watson  approximation  smoothly;  whereas, 
one  can  not  obtain  the  Debye  approximation  from  the  Watson  approxi¬ 
mation  for  small  values  of  a  |$|  which  is  the  region  where  Debye 
approximation  Is  valid.  Hence,  one  concludes  that  the  represen¬ 
tation  given  in  part  (la)  is  superior  to  the  one  given  In  part  (lb) 
of  this  appendix. 

c.  Creeping  Wave  Formulation 

Even  though  one  can  directly  obtain  the  creeping  wave  modal 
representation  for  (B-8),  It  Is  more  convenient  for  later  purposes 
to  use  (B-7)  with  p*p'*a  and  $'*0  (i.e.,  both  the  source  point  and 
the  observation  point  are  located  on  the  surface  of  the  Impedance 
cylinder  and  arc  length  Is  measured  from  the  source  point).  This 
specialized  form  of  (B-7)  Is  given  as  follows: 
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G"  *- 


k  Th  d| 

— je 


U) 


H'z'(ka) 


>|p| 


(B-34) 


i(2)  * 


Oca) 


+  3H(2) 

K  V 


(ka) 


As  mentioned  earlier,  the  creeping  wave  modal  representation 
is  valid  for  large  values  of  arc  length,  a|$|  and  it  is  obtained 
by  evaluating  (B-34)  via  the  calculus  of  residues.  To  this  end, 
one  would  use  the  Watson  approximation  for  the  Hankel  functions 
in  (B-34),  leading  to 


h  _  e-Jka|<H 
4  nm 


-®-je 


do 


W2(a)  e-jkam|*| 

w2(°)-ahmk"  ' -2(°) 


(B-35) 


where  the  use  of  (B-23)  Is  made  In  the  process  going  from  (B-34)  to 
(B-35).  The.locatlon  of  the  poles  of  the  Integrand  In  (B-35)  will  be 
denoted  by  oJJ  and  they  are  the  solutions  of 

w^(o)  -  a^mk^Wpfa)  =  0.  (B-36) 


The  solution  of  (B-36)  and  the  extensive  numerical  tables  of  the 
values  of  a"  for  different  ranges  of  the  parameter  (a.m)  are  given 
elsewhere J20],  and  they  will  not  be  Included  here  for  the  sake  of 
brevity.  From  the  calculus  of  residues,  it  can  be  shown  that 


Gh 


--2  id 


e-jka |$| 
4-imi 


I  e 

Pa0  Op+(ahmk"^)2 


(B-37) 


where  (B-36)  and  the  differential  equation  for  the  Airy  function 
(l.e.,  wi(oj-  w2(a)*0)  Is.used  In  the  process  of  obtaining  (B-37). 
One  notes  that  the  root  a"  (l.e.,  p=0  case)  Is  chosen  such  that 
It  corresponds  tohthe  surface  wave  for  the  Inductive  surface;  where¬ 
as,  there  Is  no  o"  root  for  the  capacitive  surface  and  the  summation 
In  (B-37),  for  this  later  case,  starts  from  p*l  (rather  than  p=0). 

II.  TMZ  or  Magnetic  Line  Dipole  Source  Excitation  Case 


_  The  magnetic  source  density  for  the  TM  case  Is  denoted  by 
M^,  and  It  Is  given  by  z 
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Md  =$Md  %  P~a)  ^ 


(B-38) 


where  M.  Is  the  strength  of  the  source.  The^tangential  magnetic 
field  generated  by  such  a  source  has  only  a  $-cpmponent  and  it  will 
be  denoted  by  H  .  H  can  be  obtained  from  the  z  directed  electric 

field  E*  generated  by  an  electric  line  source  T  =  z  I  p(P“a^P-^ 

where  I  Is  the  strenth  of  the  line  source.  E*  satisfies  (B-2)  and 
(B-3)  with  the  exception  that  Y  M  is  replacedjby  Z  I  and  is  sub¬ 
stituted  by  cu  where  a  is  given  in  (A-31).  E*  also  satisfies  the 
radiation  condition  for  an  eJul  time  dependence. 

With  above  changes,  (B-4)  becomes 


E^p',*’)  =  -jkZQI  6S(P'  | p •  4>) 


(B-39) 


Using  the  procedure  of  [26]  ,  one  can  show  that 


s5(p',*'iP,»)  .  -  i  r  *  *> 

-«-jc 


/ 1  \  Q_H^(ka) 

Hy1  (kP' )  -  — - H(2)(kp') 

QsHv2)(ka) 


•  H|)2^(kp)e’Jv^”^'  I  ;p>p' (B-40) 

where  Qs  ■  afiby  *  IT  • 

Actually  the  term  cosv( n- ( <|>-d> '  | )/sinvn  of  (B-6)  is  present 
in  the  construction  of  the  complete  cylinder  Green's  function  for 

the  TM  case.  However,  only  the  term  j  ^  on  the  right 

hand  side  of  (B-6)  Is  of  interest  In  this  analysis;  this  is  the 
only  term  which  is  included  in  (B-40)  just  as  in  (B-7).  The  case 
for  p'>p  can  be  obtaiogd  from  (B-40)  by  interchanging  p  and  0^. 

The  magnetic  field,  generated  by  I  can  be  found  from  VxE*z=- 
-juy  R  .  The  electric  field,  E  z  generated  by  Rd  of  (A-38)  can 
now.Be  found  from  reciprocity  tneorem  via  the  equation  Iz*E=-M. 

$*Tr  as  follows 

eJ(p.*|p\0)  *  -MdasGe(p,4>|p',0)  P>P'  (B-41a) 

EJ(P^IP'.O)  *  Mdk  Ge(p,4>|p',0)  p'>p  (B-41b) 
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where 


Ge(p.4>|p'.0)  r  dv[H(l),(kp')  - 

s  -®- je  L  *  i 

H<2),<kp)] H^IKp)  .-^lOl  i 


QsH(vU(ka) 
QsH^}(ka) 


and 


P  >P 


Q.H^dca) 


fie(p^|p'.0)  *  J  £  dJ^kp)  --5-JU - 

8  —  it  Lv  QSHv  (ka) 


H!,2)(kp)  H(12),(kp')e'JvM  ; 


P  >P 


(B-42a) 


(B-42b) 


Equations  (B-42a)  and  (B-42b)  can  further  be  simplified  for  p=a 
and  p'*a,  with  the  aid  of  the  following  wronskln  relationship 


H(v2)(ka)H(v1),(ka)  -  ' (ka)H(v1} (ka)  =  ijL  , 

yielding: 


se(qV)  -  S.  <H 


-oo-je 


.(2), 


H^'(ka) 


QsH;c'(ka) 


Q+  *  1 1m(p=a+e ;♦) 

e+0 


(B— 43) 


(B-44a) 


t)  *  (p*a;$) 

Q‘  *  11m  (p'*a+e;0) 

e+0 


(B-44b) 
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#  . 


,k. 


Making  use  of  the  other  Maxwell's  curl  equation,  H^=-(jajy  )"* 
[vx"E+WJ  ,  It  can  be  shown  that  0 

H  (p,«)  *  -JkYoMdGs(p,4>|a,0)  (B-45) 

where 


Gs(p ,♦! a,0) 


^s  i  7 
T  SI  J 


Hl2)(ka)+r  Ht2)<ka> 


e-J  V|4>|;  p>3  (B.46) 


the  prime  over  the  Hankel  function  on  (B-46)  indicates  derivative 
with  respect  to  the  argument.  Equation  (B-46)  is  specialized  to 
the  case  where  the  source  point  is  on  the  surface  of  the  cylinder 
such  that  p'=a  and  4>*  =0  and  the  observation  point  is  exterior  of 
the  cylinder  (i.e.,  p>p'=a). 

The  expression  in  (B-46)  can  be  represented  as  a  perturbation 
of  the  fields  on  the  perfectly  conducting  cylindrical  surface  as 
follows: 

GS  =  G*  +  Gp  (B-47) 


where  G*  is  the  Green's  function  for  the  magnetic  line  dipole  source 
of  (B-3§)  when  this  source  is  placed  on  a  perfectly  conducting  cylin 
der.  G*  is  obtained  from  (B-46)  upon  taking  the  limit  as  a  tends 
to  infinity.  Thus,  s 


G 


s  _  1 
C  Tit 


H^'fka) 

H^2)(ka) 


e~j  I 


(B-48 ) 


Following  an  analysis  which  is  based  on  some  earlier  work  by 
Hassejian  and  Ishimaru [19],  and  Pathak  and  Huang  [2]  one  obtains 
the  following  approximation  for  (B-48): 


H{2)(ka|$| 

kaT^l 


(B-49) 
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Gp  is  the  perturbation  Green's  function  and  it  is  given  as 
follows: 


Pv(kajots)  [pv(ka;0)]’2 


(B-50) 


where  as  before 


Pv(Ka;as) 


H^2)(ka) _ 

H(Z)'(ka)+^  H<2>(ka) 


(B-51) 


One  notes  that  in  both  (B-48)  and  (B-50),  the  observation  point, 
as  well  as  the  source  point,  is  brought  on  to  the  surface  of  the  cy¬ 
linder  (i.e.,  p'=p=a). 

As  in  the  TE  case  presented  in  section  I  of  this  appendix, 
the  perturbation  Green's  function,  Gs  of  (B-50)  can  also  be  evalu¬ 
ated  asymptotically  for  the  different  ranges  of  the  parameters  in¬ 
volved  via  the  Debye,  Watson  and  creeping  wave  representations  for 
small,  moderately  small  and  large  values  of  the  arc  length,  ka|$|, 
respectively.  In  part  ( 1 1  a)  of  this  appendix,  the  Debye  approxi¬ 
mation  of  (B-50)  is  examined.  The  Watson  approximation  and  the 
creeping  wave  representation  of  (B-50)  is  discussed  subsequently 
in  parts  (lib)  and  (lie),  respectively. 

a.  Debye  Approximation 

The  series  expansion  of  P  (ka;ctj  resulting  from  the  Debye 
approximation  can  simply  be  obtainedsfrom  (B-15b)  by  changing  Qh 
to  <£,  which  yields  n 

Pv(ka;a$) - 1 -  1  +  ->L_  - 1 - + 

Slnn+Ja.lc'1  1  2sin  "  s1n"+*>sk' 

(B-52a) 


It  can  also  be  shown  that  the  term  Py2(ka;0)  in  (B-50)  has  the  fol¬ 
lowing  form  when  it  is  expanded  via  the  Debye  approximation: 


P~2(ka;0)  ~  -sin^ 


.  J_ 

ka 


sin^n 


where  as  before  cosn^»/ka. 


(B-52b) 
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Substitution  of  (B-52)  into  (B-50)  and  keeping  the  terms  only 
up  to  the  order  of  1/v  ,  leads  to 


(B-53) 


(B-54) 


k 


;  d5(k^- 


-oo-je 


c2) 


J??  +  2jas 


e-J'Ca|«l>| 


(B-55) 


A  change  of  variable  (i.e.,  ^  )  has  been  introduced  in  (B- 

f~2  p 

51)  and  it  is  recognized  that  cosv  =  v/ka  and  sinri=Jk  -5  /k.  Com¬ 
parison  of  (B-54)  and  (A-38)  reveals  that  these  two  equations  are 
the  same  as  long  as  one  uses  the  arclength,  a]<j>|  as  the  distance 
between  the  source  and  the  observation  point.  This  distance  was 
denoted  by  | x|  in  (A-38).  Hence  the  same  argument  presented  in 
Appendix  I  can  also  be  used  here  for  the  evaluation  of  (A-38),  this 
results  in 


(B-56) 


where,  as  in  Appendix  I, 
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(B-56a) 


ss‘rJi!  “I  ■  -\s\ 


(B-56b) 


In  order  to  evaluate  Gp  ,  one  rewrites  (B-55)  as  follows: 

=  k  A  -  ~  4— ^ /"iE  d5  ,e~— - -  (B-57) 

Pi  Llc 

Comparison  of  (B-17)  and  (B-57)  reveals  that  the  integrals  involved 
are  the  same  as  long  as  ot  in  (B-17)  is  replaced  by  a  .  Hence  (B- 
57)  becomes 

„s  .1  A  d 

-  ?J  I1’  “s  IS-Jl— 


(B-58) 


where  the  notation,  ah-*-ac  in  G^(ah+oi  )  of  (B-58)  indicates  that  ah  in 
G"  of  (B-19)  will  bereplacedpbyna  .  Or,  more  explicitly,  (B-58) 


can  be  rewritten  as 


d  \f  k  te'Jkal»l 


sBvr  ^ 


_ _  -kaUlG: 

Q(j5Jka|4>|)  e  ;  a  So 


(B-59) 


The  derivative  term  in  (B-58)  or  in  (B-59)  can  be  evaluated  with  the 
aid  of  (B-22),  resulting  with 

kG^  (at+a.)1  .2  ,i,i  *  - 

-  d  po  h  s  ...  ,,  k“s  -J8sal*l  J“s  Jkald.1 

“scST  -E - (0Bsa|*|  +  1)  -g-  e  £s  +  F”  JTT" 


•  e-Jkal*l  *  e-Jkal*l  Q(jBsJi^|)e"k^a|<H:  ^0 

kJiT 

(B-60) 
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b.  Watson  Approximation 


For  moderately  small  distances  from  the  source,  one  again  em¬ 
ploys  the  Watson  approximation.  In  order  to  apply  this  approxi¬ 
mation,  one  again  starts  by  approximating  the  Hankel  functions  by 
the  Fock-type  Airy  functions  as  in  (B-23).  Following  the  arguments 
presented  in  part  (lb),  the  Watson  approximation  for  P  (ka;  a  ) 
of  (B-51)  can  be  readily  written  from  (B-26)  upon  changing  a,sto 
a  ;  that  is  h 


Pv(ka;cts) 


1 

ka 


+ 


(B-61) 


-1  3 

where  as  before  o=m  (v-ka)  and  2m  =ka.  It  can  also  be  simply  shown 
that  P'^kajO)  has  the  following  form 


P~2(ka;0)  _  (m^Vo)2 


i  i  .  n/i _ : 

(mlJB)3  H  (  [ka]2. 


(B-52) 


Substitution  of  (B-61)  and  (B-52)  into  (B-50)  yields 


Gs  „  Gs 
P  P, 


+  GS  +  0  f - - o 

pi ka  v  fka]2 


where 


'S  _  e 


-  jka  1 4>  | 


pQ  2k  4  k  2ttj 


(I  L_  J; 

Jk  2id  J 


-j«  /s+jjj7ka 


and 


(B-63) 


isa|<j>! 

ds  ------  e2  (B-64) 


Gs  =  1  e-jkaM  L_  d  Js  +  Pj/j7kas  2  sa\$\ 


(B-65) 


In  the  process  of  obtaining  (B-64)  and  (B-65),  one  employs 
a  change  of  variable  s  =-j2^  -kj;  and  naturally  m ~^Ja=  Jjs/k. 
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Evaluation  of  (B-64)  via  the  method  of  inverse  Laplace  trans¬ 
formation  yields 


Comparison  of  (B-67)  with  the  results  obtained  from  the  Debye 
approximation  for  Gp  given  in  (B-56)  reveals  after  some  manipu¬ 
lation  that  (B-56)  reduces  to  (B-67)  for  large  values  of  the  arc- 
length,  ka|<j>  and  small  values  of  a_;  whereas,  one  needs  to  conjecture 
that  s 


should  be  replaced  by 

(  i  H<2)(ka|*|)'| 

\  ka | <j> |  ) 

in  order  to  obtain  (B-56)  from  (B-67).  (It  is  noted  that  the  pre¬ 
vious  expression  is  the  large  argument  approximation  of  the  latter 
one. ) 

Hence,  one  concludes  that  (B-56)  is  a  more  general  result  than 
(B-67)  because  (B-56)  goes  to  (B-67)  for  large  values  of  the  arc- 
length  where  (B-67)  is  valid;  whereas,  (B-67)  does  not  readily  go 
to  (B-56)  for  small  values  of  the  arclength  where  (B-56)  is  valid. 

In  order  to  evaluate  G®  of  (B-65),  one  rewrites  the  equation 
as  follows  P1 


7saM 

e _ 

ysps+jJiTk  eg. 


(B-68) 


The  integrals  in  (B-28)  and  (B-68)  are  the  same  with  the  exception 
that  in  the  former  is  replaced  by  a$  in  the  latter  integral. 


Hence,  (B-68)  becomes: 


Pi 


1  / 

\  d  ^ 

kGS0(V“s> 

5  *\J 

l  K 

(B-69) 


where  the  rotation  cth-*-a$  in  (B-69)  indicates  that  ot^  in  Gp  of  (B- 
30b)  will  be  replaced  by  a$  in  (B-69)  as  before.  Or  more  Explicitly 


e-Jka|<|>| 

"  J n 


— -kal*|b! 
Q(jbs/ka|$|)e  s 


a_<0 

s 


(B-70) 


Comparison  of  (B-59)  and  (B-70)  reveals  that  the  latter  one  is 
the  special  case  of  (B-58)  for  small  values  of  a  One  notes  that 
the  derivative  term  in  (B-69)  or  (B-70)  can  be  obtained  from  (B- 
33)  by  changing  to  a$. 

c.  Creeping  Wave  Formulation 


With  both  the  source  and  the  observation  points  on  the  surface 
of  the  cylinder  (i.e.,  p=p'=0),  (B-46)  becomes 


rs  .  Jas  1 


Hr  C  -fe) 


H 


(2)* 


(ka) 


Hl2)(ka)+r  C'(ka) 


cl 


(2), 


5-Jv|4>|(8_71) 


A  creeping  wave  or  residue  series  representation  will  now  be 
given  for  (B-77).  In  the  process  of  obtaining  a  residue  series 
expansion,  one  uses  the  Watson  approximation  of  (B-23)  into  (B-71), 
yielding 
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The  poles  of  the  Integrand  in  (B-72)  will  be  denoted  by  and  they 
are  the  roots  of  following  characteristics  equation  p 


w£(op)  -  asmk-1w2(op)  =  0 


(B— 73) 


As  indicated  earlier  for  cr^,  extensive  numerical  tables  for  the 
values  of  o*  for  differentpranges  of  the  parameters  (a  m)  are  also 
given  by  Logan  [2CQ  ,  and  they  will  not  be  included  heresfor  the  sake 
of  brevity.  From  the  calculus  of  residues  and  the  differential 
equation  for  the  Airy  functions  (i.e.,  wJJ(a)-  w?(cr)=0),  one  can 
show  that 
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In  the  TM  case,  root  (I.e.,  p=0  case)  corresponds  to  the  surface 
(or  Elliot)  wave  for  the  capacitive  surface,  whereas,  there  is  no 
cr  root  for  the  inductive  surface  and  the  summation  in  (B-74),  for 
this  latter  case,  starts  from  p=l  (instead  of  p=0). 
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